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The electric quadrupole transitions between 0"^, 2’*', and 4'*' states in are investigated in a 
3a model. The three-body wave functions are obtained by means of the hyperspherical adiabatic 
expansion method, and the continuum is discretized by imposing a box boundary condition. Cor¬ 
responding expressions for the continuum three-body (3a) bremsstrahlung and photon dissociation 
cross sections are derived and computed for two different a — a potentials. The available exper¬ 
imental energy dependence is reproduced and a series of other cross sections are predicted. The 
transition strengths are defined and derived from the cross sections, and compared to schematic 
rotational model predictions. The computed properties of the resonances suggest that the two 
lowest bands are made, respectively, by the states {O)*", 2)'", 4^} and {0^, 2^, d)*"}. The transitions 
between the states in the first band are consistent with the rotational pattern corresponding to 
three alphas in an equal sided triangular structure. For the second band, the transitions are also 
consistent with a rotational pattern, but with the three alphas in an aligned distribution. 

PACS numbers: 23.20.-g, 21.60.Gx, 21.45.-v 


I. INTRODUCTION 

The structure of the spectrum has attracted a lot 
of attention along the years. In this nucleus only two 
bound states exist, the 0)*" ground state and the first ex¬ 
cited state, with angular momentum and parity 2)'’ and 
excitation energy of 4.44 MeV. We use Jf to represent 
the state with angular momentum J and parity tt. 
Already in the 50’s, F. Hoyle predicted the existence of 
a Oj resonance at an excitation energy of 7.65 MeV, as 
a requirement in order to explain the known abundance 
of carbon in the Universe [I|. The existence of such res¬ 
onance was experimentally confirmed just a few years 
later Q] . Its properties have recently been thoroughly re¬ 
viewed Q. Since then, many other resonant states have 
been observed in (dl. iQ-HOjl . Among them, one of the 
most elusive ones has been the second 2^ state, which 
was predicted by different theoretical methods to have 
an excitation energy of around 10 MeV EMI. Recent 
experiments have confirmed the existence of this 2^ res¬ 
onance at the expected energy Biii- 

Together with the O'*' and 2+ states mentioned above, 
a resonance is known to be at an excitation energy of 
14.1 MeV B0- Other 4+ states have also been obtained 
numerically. In fact, the known resonance is usually 
found to be the 4^ state, and a first resonance at a lower 
energy is often obtained (for instance, at 11 MeV in E3| 
or 10.5 MeV in EH). Experimental evidence of the 4^ 
state was reported in Ref.E^, and its excitation energy 
was given to be 13.3 ± 0.2 MeV, only about 1 MeV lower 
than the 4^ state. 

The appearance of the {0+,2+,4+} sequences in the 
energy spectrum suggests the rotational character of 
these states. Therefore at least two rotational bands 
could exist in one of them sitting on the ground 
state, and another one sitting on the Hoyle state. In fact, 
it is becoming very common in the literature to refer to 


these states as rotational states El El- 

The rotational sequences of states are not necessarily 
{O'*", 2+, 4+}, since the underlying intrinsic shapes differ¬ 
ent from axial and R 2 {t^) symmetry may provide addi¬ 
tional states in a given band. This is the classical knowl¬ 
edge that the quantum numbers specifying the states in 
a rotational band directly carry information about the 
symmetry of the intrinsic state. The details of conse¬ 
quences of the D^h symmetry (equilateral triangle) for 
the 3 q: spectrum of was formulated and discussed 
in general in El- Evidence for that symmetry was pro¬ 
vided in El by comparing energy sequences and tran¬ 
sition probabilities. We emphasize that also this model 
treats the resonances as bound states where all strengths 
are collected in the corresponding bound-state wave func¬ 
tion. 

A simple rotational sequence is also observed in ®Be, 
where the ground state and the two first excited states 
follow as well the angular momentum and parity sequence 
0+, 2+, and 4+ (also 6'*' and S’*' states with large widths, 
are found numerically as poles of the 5-matrix). All the 
states in the ®Be spectrum are unbound (resonances), 
with the ground state only about 0.1 MeV above the 
threshold for emission of two alpha particles. There¬ 
fore, the question arises about a possible rotational (or in 
general collective) character for a sequence of continuum 
states of considerable width. T his p roblem was recently 
investigated in previous works EMt 21|. This was done 
by computation of the electric quadrupole cross sections 
for bremsstrahlung emission after transitions between the 
®Be-states, which were described as two-body systems 
made of two alpha particles. These calculations required 
a careful treatment of the continuum wave functions and 
clarification of the definition of the cross section for tran¬ 
sitions between states with a none well-defined energy. 

From these ®Be cross sections it is possible to extract 
the transition strengths, which for the case of transitions 
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within a rotational band must follow well-established 
rules based on the assumption that all the states have the 
same intrinsic spatial structure (rigid rotor). The main 
conclusion was that the computed transition strengths do 
not behave as expected for states in a rotational band. In 
fact, when increasing the angular momentum of the ini¬ 
tial states the transition strength was found to decrease, 
which is precisely the opposite to the prediction of the ro¬ 
tational model [^. Nevertheless, allowing the separation 
between the two alphas in the ®Be-resonances to change 
according to the computed root mean square radius, a 
very nice agreement between the computed strengths and 
the rotational model predictions was observed. 

The purpose of the present work is to extend the ®Be 
continuum investigations to the spectrum of and 
check if the {0)*', 2('', d)*"} and {Oj, 2^,42} sequences of 
states follow the behavior predicted by the rotational 
model. For this aim we shall compute the electric 
quadrupole q-emission bremsstrahlung cross sections for 
the different transitions between the states. These cross 
sections and the transition strengths are not well-defined 
quantities when involving resonances with finite width. 
Suitable definitions are necessary to be specific. We 
shall work entirely within the strict two- and three-body 
framework. This means that the constituent particles are 
inert a-particles, and the inter-particle nucleonic Pauli- 
principle will be accounted for without use of the intrin¬ 
sic microscopic structure. This model, used throughout 
the present paper, can naturally be referred to as an a- 
particle model, in contrast to microscopic a-cluster mod¬ 
els where the intrinsic nucleonic composition of the a- 
particle is the basic ingredient [l^ [Ij . This procedure 
follows very closely the one described and tested in [ 2 l| 
for the two-body system of ®Be. In this way we can test 
the validity of the generalization to three-body systems. 

The structure of will be approximated as a three- 
alpha system and calculated by use of the hyperspherical 
adiabatic expansion method [22|. The continuum spec¬ 
trum will be discretized by imposing a box boundary con¬ 
dition in the hyperradial coordinate, and the correspond¬ 
ing wave functions will be computed on the real energy 
axis, without any particular treatment of the resonances. 
All states are equally treated as continuum states charac¬ 
terized by energy, angular momentum, and parity. A res¬ 
onance would only be a distribution of these discretized 
continuum states. 

Two different alpha-alpha potentials will be used in 
our calculations, the Ali-Bodmer and the Buck poten¬ 
tials [mil. These potentials are parametrized to re¬ 
produce a —a scattering properties and consequently the 
two-body wave functions must be the same at large dis¬ 
tances. However, both small distances and high partial 
wave properties can differ drastically. All the available 
versions of the Ali-Bodmer potential reproduce the ex¬ 
perimental s- and d-wave phase shifts, but only versions 
’d’ and ’e’ reproduce also the ones for £ = 4 [2^. For 
this reason we have chosen to use the version “d” of the 
Ali-Bodmer potential which has no bound states in con¬ 


trast to the Buck potential [2^. Nevertheless, as shown 
in Ref. [ 2 ^, Ali-Bodmer and Buck potentials give rise to 
very similar phase shifts for i = 0 , 2 ,4, 6 , and 8 . 

The main difference between the two potentials is that 
the Buck potential contains two s-wave and one d-wave 
Pauli-forbidden bound states. This fact implies a very 
different short-distance behaviour for the 0 "^ and 2 + wave 
functions in ®Be, since the corresponding radial wave 
functions show a different number of nodes depending 
on the potential used. For this reason, observable quan¬ 
tities sensitive to the short-distance structure, like the 
photon emission, were expected to provide information 
about the underlying two-body potentials. As shown in 
Ref, . to our surprise, we found no significant differ¬ 
ences in the computed bremsstrahlung cross sections for 
®Be between these potentials. 

At the three-body level things are very different. For 
instance, the ground state in is always the lowest 
computed O"*" state, no matter which of the two potentials 
is used. This is because when using the Buck potential we 
have two options. Either we construct the phase equiv¬ 
alent alpha-alpha potentials [ 2 ^ or we exclude the Pauli 
forbidden adiabatic potentials before computing the ra¬ 
dial three-body wave functions [1^. In both cases the 
effective two-body potential shows a short-distance core 
repulsion, similar but not identical to the one in the Ali- 
Bodmer potential. As a consequence, the structure of 
the three-body system looks very much the same in all 
the cases. Therefore, we could off hand expect for 
a dependence on the potential even smaller than in the 
®Be case. This point will be tested in our calculations. 

The paper is organized as follows. In Section II we 
describe the theoretical background needed to compute 
the cross sections. The particularization of the general 
expressions to the three-body case is given in Section 
III, which is divided into three subsections devoted to 
the three-body wave functions, the transition matrix el¬ 
ement, and the transition strengths, respectively. In sec¬ 
tion IV we focus on the three-alpha system, describing 
the details of the two-body potentials used in the cal¬ 
culation and discussing the corresponding computed 
spectra. The electric quadrupole 7 -emission cross sec¬ 
tions are described in Section V. Finally, in Section VI 
we give the computed A 2 -transition strengths and com¬ 
pare with the predictions from the rotational model. We 
finish the paper with the summary and the conclusions. 
Some derivations and expressions not essential for the un¬ 
derstanding of the paper, but important in order to make 
the paper self-contained, have been collected in three ap¬ 
pendices. 


II. CROSS SECTION EXPRESSIONS 

The photo-dissociation cross section for the breakup 
of a bound system A with angular momentum Ja into a 
three-body system with angular momentum J and three- 


3 


body energy E (A + 7 —>■ a + 6 + c) is given by 




(27r)^(A + l) 
A[(2A + 1)!!]2 



dE 


{Ja —>■ J), 


( 1 ) 

where A is the multipolarity of the electromagnetic tran¬ 
sition, E^ is the photon energy (E^ = E + \Ba\, where 
Ba is the binding energy of A), and 


dE 


{Ja J) 


1 

2t/yi -|- 1 




^ 5{E-Ei), 


( 2 ) 

where we have assumed that the continuum spectrum de¬ 
scribing the final three-body system has been discretized 
and the index i runs over the discrete continuum states. 

The wave function of the discrete state i, with energy 
(i) 

Ei, is denoted by $y, and is the wave function of 
the bound state A. The electromagnetic operator with 
multipolarity A is denoted by 0\, which for the case of 
electric transitions reads: 


dx^=eY,Z^r^Yxi^{n,), (3) 

i 


where i runs over the charged particles in the system, 
each of them with charge eZi, and is the center of 
mass coordinate of particle i, whose direction is given by 
the angles Vli. 

As described in appendix m the photo-dissociation 
cross section a{^\E) corresponding to the A -|- 7 —>■ 
a-|- 6 -|-c process, and the one corresponding to the inverse 
reaction a + h + c^ A-I -7 (denoted now as (A)), are 
related through the expression: 

a^^\E) _ , 2{2Ja + 1) ?,2tt fE^V 

“''■(2Ja + l)(2Jb + l)(2J, + l)^ V'fcj ’ 

(4) 

where v is the number of identical particles in the three- 
body system, Ja, Jb, and Jc are the total angular mo¬ 
menta of particles a, b, and c, respectively, and k is 
the three-body momentum, which is defined as k = 
yj2mEjW. The mass m is the normalization mass used 
to define the Jacobi coordinates, which are the coordi¬ 
nates usually employed to describe the three-body sys¬ 
tem [^ . 

As seen in Eq.Q, the cross section a^^\E) corre¬ 
sponding to the radiative capture reaction a + b + c ^ 
A J- 7 has dimensions of length to the fifth power, which 
corresponds to a surface in the six-dimensional space re¬ 
quired to describe the incoming three-body system. This 
cross section depends through k on the normalization 
mass m. This is related to the fact that when using the 
hyperspherical coordinates (which are constructed from 
the Jacobi coordinates), the radial coordinate, the hyper¬ 
radius, depends also on m. Therefore, a given value of 
the hyperradius will correspond to different three-body 
geometries (different relative distances between the three 
particles) for different choices of m. As a consequence, 


the flux of incoming particles through a given hypersur¬ 
face of radius p will depend on m, and therefore also the 
cross section, which is the outgoing flux of particles nor¬ 
malized with the incoming flux (see appendix . Note 
that the well-defined physical observable, independent of 
the choice made for the normalization mass, is the reac¬ 
tion rate (see appendix iBl). 

If we now replace the bound state. A, by a contin¬ 
uum state with energy E' and momentum J' , the photo¬ 
dissociation cross section for the process A-I -7 —>• a-|- 6 -|-c 
given in Eq. o can be easily generalized as described in 

Ref. [HI: 


da. 


(A) 


dE' 


■{E) = 


(27r)3(A-bl) fE. 


A[(2A-bl)!!]2 \hc 


2A-1 


dEdE' 


where Ej = E — E' is the photon energy and 


{J' ^ J), 

(5) 


dEdE’ 

1 

2J' + 1 


(J' ^ J) = 


( 6 ) 


^ I|dy !<!>«) d{E - E,)6{E' - A'), 




where the initial continuum states have also been dis¬ 
cretized, and the index j runs over the initial continuum 
discrete states with energy A' and wave function It 
is important to note that the summation over i and j in 
the equation above is not unrestricted, but limited to the 
initial and final energy ranges of experimental interest. 
We will come back to this point later on. 

Einally, according to Eq.®, we have that the cross 
section for the continuum to continuum reaction a + b + 
c —>■ A J- 7 , with initial and final energies E and E', and 
initial and final angular momenta J and J' is given by: 


da^^'f 

dE' 


{E) = 


(7) 


2{2Ja + 1) 


’{2Ja + l){2Jb + l){2Ja + l) At® 


327r / E.y 
he 


’ <i7'’ 




As discussed in Ref. [2l| , the total bremsstrahlung cross 
section, as a function of the incident energy E, is obtained 
after integration over E' (or over the photon energy A-,,). 
In our description using discrete initial and final contin¬ 
uum states, the integral over E' is trivially made thanks 
to the 5{E' — E'j) function in Eq.®, and just a summa¬ 
tion over the index j running over the final discrete states 
remains: 


r(A) 


{E) = 


{2Ja + l){2Jb + l){2Ja + l) \hcj 


E,\ 


2\+l 

( 8 ) 


16 1 (27r)5(A-H) 
TT AC® A[(2A + 1)!!]2 


^ (<i>«|| 6 v||$y?) 6{E-E,). 




As anticipated below Eq. ® and stated in Ref. [T^ , the 
computed cross sections should be obtained in close anal¬ 
ogy to the experimental setup, where only a finite range 
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of final relative energies is measured, usually around a 
resonance in the final system. This means that the in¬ 
tegral over E' has to be performed only over this pre¬ 
cise energy range. In our language of discrete continuum 
states, this means that the summation over j in the equa¬ 
tion above runs only over the final discrete states whose 
energy E'^ is contained in the chosen final energy window. 
It is obvious then that the total bremsstrahlung cross sec¬ 
tion depends on such final energy window, see [20l l2l| . In 
other words, definitive statements about resonance prop¬ 
erties must take into account that these continuum states 
have an energy at best defined with an accuracy of less 
than its width. 


III. THREE-BODY INGREDIENTS 

In order to compute the cross section given in Eq.([5|) 
we first need the wave functions for the three interacting 
particles in states with given energy and specified an¬ 
gular momentum and parity. Then we need the transi¬ 
tion probability from one state to another. This requires 
specific matrix elements which by combination of vari¬ 
ous kinematic factors provide cross sections and through 
proper definitions also the related transition strengths. 

A. Three-body wave functions 

In this work we shall construct the three-body wave 
functions using the hyperspherical adiabatic expansion 
method described in Ref.[^. In this method the wave 
function is expanded in terms of a complete set of angular 
functions where J is the total angular momentum 

of the three-body system: 

= fn {p)€{p. ( 9 ) 

P TJ. 


where p = a = arctan(x/?/), and 

are the angles defining the directions of x and y, which 
are the Jacobi coordinates used to describe the system. 
Writing the Schrddinger equation in terms of these co¬ 
ordinates, they can be separated into angular and radial 
parts: 


In Eg. (ITT]) E is the three-body energy, and the coupling 
functions Pnn' and Qnn' are given for instance in p^ . 
The potential Vstip) is used for fine tuning to take into 
account all those effects that go beyond the two-body 
interactions. 

It is important to note that the angular functions used 
in the expansion in Eq.(|ni) are precisely the eigenfunc¬ 
tions of the angular part of the Schrddinger (or Ead- 
deev) equation(s). Each of them is in practice obtained 
by expansion in terms of the hyperspherical harmonics, 
see Ea. (IC2l) . Obviously this infinite expansion has to be 
truncated at some point, maintaining only the contribut¬ 
ing components. 

The eigenvalues A„(p) in Ea. dTUI) enter in the radial 
equations (ED as a basic ingredient in the effective radial 
potentials. Accurate calculation of the A-eigenvalues re¬ 
quires, for each particular component, a sufficiently large 
number of hyperspherical harmonics. In other words, the 
maximum value of the hypermomentum, Kmax, for each 
component must be large enough to assure convergence 
of the A-functions in the region of p-values where the 
fnip) wave functions are relevant for the calculation of 
the electromagnetic operator matrix element. 

Finally, the last convergence to take into account is 
the one corresponding to the expansion in Eq. (P|) . Typ¬ 
ically, for bound states, this expansion converges rather 
fast, and usually three or four adiabatic terms are already 
sufficient. 

In our calculations the continuum is discretized by use 
of a box boundary condition. This means that the radial 
wave functions fn{p) are imposed to be zero at a given 
maximum value of p, which is typically taken equal to 
a few hundreds of fm {pmax = 200 fm in our calcula¬ 
tions). No distinction is made between resonances and 
ordinary continuum states. Therefore, in order to place 
the initial and final energy windows matching the res¬ 
onance energies, it will be necessary to have some pre¬ 
vious information about what these energies are. It also 
implies that the computed result simultaneously includes 
on-resonance as well as continuum background contribu¬ 
tions. 


B. Transition matrix element 




A>; 


2mp^ 

~1^ 


{Vju+V.k + V.j)(t>i = Xu{p)cl>l ( 10 ) 


and 


Once the initial and final three-body wave functions 
are computed, we can now obtain the square of the tran¬ 
sition matrix element 


(P 2to 1 \ 15A 

~d^ + -?)^(^36(p) -E) + ^ i^Xnip) + ^ j 

+ J2n' “ Qnn'^j fn'ip) = 0, 


fnip) 

( 11 ) 


where Vjk, Vik, and Vij are the two-body interactions 
between each pair of particles, A^ is the hyperangular 
operator, see Ref. [2^, and m is the normalization mass. 


K$«||a^ll$o;))|2 ^ (2j +1) ^ 

fiM' 

( 12 ) 

which enters in Eq.(l8]), and where Ox^ is the electric 
multipole operator given in Eq.([3]). 

From the definition of the Jacobi coordinates. Ref. [2^, 
it is not difficult to see that the vector giving the 
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position of particle p can be written as: 


I m{mq + nis) 

mp{mp + mq + nis) 


Vp 


(13) 


where rup, niq, and nis are the masses of the three parti¬ 
cles and Up is the Jacobi coordinate defined between par¬ 
ticle p and the center of mass of the other two. Therefore, 
using Eq. we get: 


p=i 


m{mq + nis) 


mp{mp 


■ rus) 


x{^f\\y^Y^{fp 




A/2 


(14) 


and the calculation of the reduced matrix element in 
Eg. dT^ requires only the calculation of the matrix ele¬ 
ment (‘hj^||yplbi(f^p)||4’j/^) for each of the three possible 
definitions of the Jacobi coordinate j/p. This reduced ma¬ 
trix element is obtained as described in appendix [Cj The 
final expression in Ea. (IC4l) is rather elaborate although 
mostly due to the Racah algebra necessary to account 
correctly for all total angular momentum couplings in 
general cases of particles with finite intrinsic spin. 

The integral over p in Eq. (|C4I) involves the radial wave 
functions f;( and f^, corresponding to the initial and 
final continuum states. Therefore the function to be in¬ 
tegrated does not fall to zero at infinity, but instead gives 
rise to an apparent divergence of the necessary integrals. 
This divergence is obviously nonphysical and mathemati¬ 
cally ill-defined until a suitable limiting procedure is cho¬ 
sen. A simple way to solve this numerical problem was 

proposed in Ref. [21| , where the integrand was multiplied 
2 2 

by the factor ^ , in such a way that the correct re¬ 
sult is obtained in the limit oi rj — 0. In practice, a value 
of 77 in the vicinity of rj — 0.01 fm“^ is enough to get a 
sufficient accuracy. 

It is important to note that the divergence mentioned 
in the previous paragraph appears due to the application 
of the long-wavelength approximation, thanks to which 
the electric field can be obtained from the charge den¬ 
sity only (Siegert theorem). This theorem leads to the 
electric transition multipole operator given in Eq.Q. In 
Ref. [ 2 ^ an extension of the Siegert theorem not relying 
on the long-wavelength approximation was proposed, in 
such a way that the divergence problem in the matrix 
element disappears. The equivalence between this proce¬ 
dure and different techniques designed to treat the diver¬ 
gence problem (the one used in this work among them) 
has been investigated in Ref. [ 2 ^. 


C. Transition strength 


the J-functions would permit to obtain the total transi¬ 
tion strength just by summation over i and j of the square 
of the reduced matrix element. This would in fact be the 
same procedure as when a bound state is involved in 
the transition. However, as discussed in Ref. [ 2 ^, an in¬ 
discriminate sum over initial and final continuum states 
makes the result rather meaningless. The information 
about resonance properties is completely washed out and, 
even worse, weighted at the wrong energies. Further¬ 
more, due to the undesired divergence produced by the 
soft-photon contribution, which appears when the en¬ 
ergy of the emitted photon approaches zero {E' —>• E 
or Ej —>■ Ei) 0, the calculation itself is pretty compli¬ 
cated. For these reasons, we shall obtain the transition 
strength as described in Ref. [ 2 ^, i.e. directly from the 
cross section in Eqs. © and ©. Two different methods 
will be used. 

In the first method the strength is obtained from the 
total (integrated) cross section. More precisely, for a 
transition between some initial and final energy windows, 
typically around some resonances in the initial and final 
states, integration of Eq. © over E and E' within those 
two windows will provide the total cross section for the 
transition. If the photon energy, Ej, were constant, this 
would immediately provide the total transition strength 
just after division by the constants multiplying the tran¬ 
sition strength. However, since this is not correct, we 
must use an average value for In particular, the 

photon energy will be taken as E^ = Ec — A' where Ec is 
the energy of the cross section peak corresponding to the 
resonance in the initial state and A' is the center of the 
final energy window (usually the resonance energy in the 
final state). Of course, this procedure assumes informa¬ 
tion about resonance positions, and it can be sensitive to 
rather small variations around a chosen E^ owin g to the 
power of (2A J- 1) for B{E\) transitions, see Ref.p^ for 
details. 

The second method exploits the fact that in the 
vicinity of a resonance the cross section of the photo¬ 
dissociation reaction A + ')^a + h-\-c takes the form 


aj{Ej) 


2 J-I-1 Trfi^c^ ri{r..y 

2{2Ja + 1) E^ {E-En)^ + ryY ^ 


where after the collision the particles a, b, and c are as¬ 
sumed to populate a resonance with angular momentum 
J, energy E/j, and width for decay into three particles 
Tfl. Ja is the angular momentum of A and T = T/j -|- T-y, 
where T.^ is the 7 -decay width of the three-body reso¬ 
nance. 

Using Eq.(|4l), we can easily obtain the expression 
equivalent to Eq. m for the inverse cross section after 
the three-body collision a + b + c ^ A-I- 7 , which can be 
written as 


In principle, calculation of the strength for the 

continuum to continuum transition, a + b-\-c^A-\-^, 
could be made directly through Eq.®, which thanks to 


a{E) = (16) 

v\[2J+l) 8(27r)2 VjiT^ 

(2Ja + l)(2Jb + l)(2Je + l) (E-E;fl)2+r2/4’ 
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where now the three colliding particles, with spins Ja, Jb 
and Jc are assumed to populate a three-body resonance 
with angular momentum J, energy and width for 
decay into three particles r^;. 

The second method uses the value of T-y in the equa¬ 
tion above in order to fit the peak in the computed cross 
section in Eq.® corresponding to the three-body reso¬ 
nance at E = Eh, in such a way that from T-,, we can 
obtain the transition strength thanks to the well-known 
expression [^ : 




87r(A -I- 1) 
A[(2A + 1)!!]2 



2A-I-1 

J a)- 


(17) 


The two methods can be compared. The first prescrip¬ 
tion depends on the windows chosen for both initial and 
final states. These choices should be made precisely to 
reproduce the conditions in a given measurement. How¬ 
ever, then the average photon energy becomes impor¬ 
tant. The second prescription relies on the behavior of 
the cross section for energies around the resonance peak. 
This dependence is assumed to have the Breit-Wigner 
shape with the advantage (and related weakness) that 
only a few points around the peak are used. The pho¬ 
ton energy does not enter as a multiplicative factor, but 
deviations from the assumed simple Lorentzian behavior 
are not accounted for. 


IV. THE THREE-ALPHA SYSTEM 

In this section we describe the calculation made to con¬ 
struct the wave functions for a system of three alpha par¬ 
ticles. We start giving the details of the two-body a-a 
potentials used in the calculation and some properties of 
the ®Be-spectrum. In the second part we summarize the 
properties of the ^^C-spectrum. 


A. ®Be properties 


We shall consider the same two a-a potentials used in 
Ref. [ 2 ^ , i.e. the Buck potential 241 and version d of the 
Ali-Bodmer potential given in Ref. [ 2 ^. The Buck po¬ 
tential has two spurious deep-lying a-a bound states for 
s-waves and one more for d-waves. These spurious states 
correspond to Pauli forbidden states. On the contrary, 
the Ali-Bodmer potential is a shallow potential not hold¬ 
ing any bound a-a state. These two potentials give rise 
to very similar ^ =0, 2, 4, 6, and 8 phase shifts. 

The spectrum of ®Be obtained with these two poten¬ 
tials is discussed in Ref. [ 2 ^. In here we only summarize 
in Table U the energies and widths of the different states. 
The first two rows show the known experimental ener¬ 
gies and widths [s^ of the O'*", 2+, and 4+ resonances 
in ®Be. The experimental values are well reproduced by 
both potentials. The computed resonances have been 
obtained as poles of the iS-matrix by use of the complex 


scaling method [3l|, The widths of the computed 
6+ and 8+ resonances are very big, comparable to their 
energies, and actually, they should not be considered as 
well-defined resonances. 

After a complex scaling calculation the complex ro¬ 
tated wave functions of the resonances fall off asymptot¬ 
ically as ordinary bound states. It is therefore possible 
to compute mean-square-radii, expectation values of r^, 
which for resonances are complex numbers, in contrast to 
the real values obtained for bound states even if the corre¬ 
sponding wave functions have been complex rotated. The 
real and imaginary parts of for each of the com¬ 

puted resonances are also given in Table [H As introduced 
by T. Berggren [H, and also discussed in Ref.js^, 
the real part of the expectation value of a given com¬ 
plex rotated operator has been attempted interpreted as 
a corresponding average value over continuum wave func¬ 
tions in a range of energies around the resonance. It is 
then tempting to associate the imaginary part of 
with an uncertainty of the resonance size arising from the 
non-zero width of the state. 


B. spectrum 

The resonances of the three-alpha system are obtained 
as described in Ref.[ll|. The method follows the hy- 
perspherical adiabatic expansion method sketched in 
Sect. IIII A[ which is used in combination with the com¬ 
plex scaling method [3l|, . The three-body resonances 
appear then as ordinary bound states with complex en¬ 
ergy, whose real and imaginary parts give the resonance 
energy and half the width of the resonance. This method 
does not make any assumption about the resonance prop¬ 
erties. For instance, the resonance decay mechanism is 
dictated by the dynamic evolution of the resonances, i.e., 
by the change in structure from small to large distances. 
In this way, the sequential and direct decay channels are 
both simultaneously taken into account, and the corre¬ 
sponding branching ratios are directly dictated by the 
resonance wave function ( 3 ^ . 

In order to reproduce the known experimental ener¬ 
gies in the ^^C-spectrum a fine tuning with a short-range 
three-body force is required. This is done by the potential 
Vsbip) introduced in the set of radial equations given in 
Eq. (fTTll . As in Ref. [ll| , we shall consider here a gaussian 
three-body potential V^b = Se~P , where the range b is 
taken equal to 6 fm, which approximately corresponds to 
the hyperradius obtained from three touching a-particles. 
This Vsb construction maintains the structure of a three- 
body state, but varying the strength, S, the energy po¬ 
sition can be adjusted to reproduce the measured value 
of the resonance. 

When using the Buck potential, due to the existence 
of Pauli forbidden two-body states, a direct three-body 
calculation gives rise to a large amount of spurious bound 
three-body states. To avoid this problem, the three-body 
calculation is made using the phase-equivalent version of 
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TABLE I. Properties of the five lowest computed resonances in ®Be. The first two rows give, when available, the corresponding 
experimental energies, Er, and widths, Tr, taken from Ref.[^. The compnted values with the Buck and Ali-Bodmer potentials 
are given by the third and fonrth rows, and by the fifth and sixth rows, respectively. All the energies and widths are given in 
MeV. The following four rows give, also for the two a — a potentials, the real and imaginary parts of v<r^>, computed with 
the complex scaling method. These values are given in fm. 


J+ 

0 + 

2 + 

4+ 

6 + 

8 + 

Er (Exp.) 

0.0918 

2.94 ±0.01 

11.35 ±0.15 

- 

- 

Tr (Exp.) 

(5.57 ±0.25)10"® 

1.51 ± 0.02 

~ 3.5 

- 

- 

Er (Buck) 

0.091 

2.88 

11.78 

33.55 

51.56 

Tr (Buck) 

3.6- 10“® 

1.24 

3.57 

37.38 

92.38 

Er (Ali-Bodmer d) 

0.092 

2.90 

11.70 

34.38 

53.65 

Tr (Ali-Bodmer d) 

3.1 ■ 10“® 

1.27 

3.07 

37.19 

93.74 

ReV< > (Buck) 

5.61 

3.51 

2.93 

2.82 

2.76 

Im%/< > (Buck) 

0.01 

1.29 

0.82 

1.44 

1.77 

ReV< > (Ali-Bodmer d) 

5.80 

3.58 

2.91 

2.70 

2.73 

ImV < > (Ali-Bodmer d) 

0.001 

1.24 

0.76 

1.40 

1.73 


the Buck potential. This potential is constructed nu¬ 
merically from the original one, and it provides a two- 
body potential with exactly the same phase shifts for 
all energies, but where the undesired forbidden bound 
states have been removed from the two-body spectrum, 
see Ref. [ 2 ^ for details. 

In Table m we give the computed energies and widths 
obtained for the 0+, 2+, and 4+ states, which are the 
states of interest for this work. The experimental data 
are taken from Refs. [1, H, d, [H, [s^ ■ The results obtained 
with the Ali-Bodmer and the Buck potentials are given. 
For each calculation the value of the strength. S', used in 
the three-body force is also given. 

For the 0+ states, if the three-body force is used to 
fit the energy of the ground state, we then get a sec¬ 
ond 0+ state (the Hoyle state) slightly bound, which 
is clearly incorrect. In order to get the correct energy 
for the Hoyle state, whose structure is expected to agree 
more with the three-alpha model than the ground state, 
it is therefore necessary to weaken the three-body attrac¬ 
tion, which makes the ground state underbound by a bit 
more than 1 MeV. When this is done, a third 0+ state 
appears at an energy of about 5.5 MeV with the Ali- 
Bodmer potential and 4 MeV with the Buck potential. 
This last energy agrees with the experimental value, al¬ 
though the large experimental width is better reproduced 
with the Ali-Bodmer potential. It is important to keep 
in mind that the resonances are computed as poles of the 
iS-matrix, and the experimental data are often obtained 
after an 77.-matrix analysis of the cross sections. These 
two different procedures can lead to sometimes very dif¬ 
ferent values for the resonance widths [s^. In Table HIl 
the two values of the strength of the three-body force S 
fitting the experimental energies of the 0)^ and Oj states, 
are given. 

For the 2+ states we proceed in a similar way. When 
the three-body force is used to fit the energy of the bound 


2 + state, we then get the second 2+ state with an energy 
and width of (1.75,0.38) MeV and (1.72,0.16) MeV with 
the Ali-Bodmer and Buck potentials, respectively. These 
values of energy and width are lower than the experi¬ 
mental values recently given in [1, |^. Again, in order to 
fit this energy we have to weaken the three-body force, 
in such a way that the 2^ resonance appears at about 
2.4 MeV. An even weaker three-body force would place 
the 2^ state at 2.8 MeV, in better agreement with the 
experimental value, but in this case the resonance would 
be too wide and more difficult to obtain through a com¬ 
plex scaling calculation (due to the large rotation angle 
required in this case). When the energy of the 2^ res¬ 
onance is placed at about 2.4 MeV, the corresponding 
width is around 1 MeV with both, the Ali-Bodmer and 
the Buck potential. This width agrees with the value 
given in Ref.@, but it is a factor of 2 smaller than the 
value given in Ref.Q. Note that the resonance widths 
given in Refs.Q and Q are obtained from the same ex¬ 
periment, although the fit of the data is obviously differ¬ 
ent. This fact emphasizes that not too much confidence 
should be placed on the comparison in Table |ll] between 
the computed and experimental widths. 

Beside the two lowest 2+ (bound and resonance) states, 
also a 2;^ state appears in the calculations. As seen in 
Table n when the energy of the 2^ resonance is moved 
from ^ 1.8 MeV to ^ 2.4 MeV, the effect on the 2^ state 
depends rather strongly on the potential used. With the 
Ali-Bodmer potential the energy moves from 3.90 MeV 
to 6.87 MeV, and the width changes quite dramatically 
from a rather small value, 0.16 MeV, to 2.01 MeV. How¬ 
ever, with the Buck potential basically only the resonance 
width changes, from 0.6 MeV to 3.6 MeV, while the en¬ 
ergy value remains quite stable, since it only changes 
from 5.04 MeV to 5.48 MeV. 

Concerning the 4+ states, the existence of a resonance 
at 6.81 MeV above the three-body threshold, with a 


















TABLE II. Calculated and measured energies Er (in MeV) and partial a-decay widths Th (in keV) of the resonances 
for different . Experimental values (labeled “exp”) are from [1, H, 0, HE, Is^]. Th e labe ls “AB” and “Buck” refer to the 
calculations obtained using the Ali-Bodmer and Buck a-a potentials specified in Sect. IIV Al The strength S (in MeV) used in 
the gaussian three-body potential Vsb is also given for each of the calculations (the strength of the three-body force is always 
taken to be 6 = 6.0 fm). The energies are measured from the three-alpha threshold. 



-S'/?, exp 

r j?,exp 

^R,AB 

Tfl.AB 

5 

Ti^fl.Buck 


S 

Ot 

_7 275(«) 

- 

-7.27 

- 

-22.6 

-7.27 

- 

-22.0 

02 "- 

0.380(“> 

0.009 ± 0.00l(“> 

0.38 

< 0.05 

-18.2 

0.38 

<0.05 

-18.0 

0^ 

4.20±0.14l*’l 

3440 ± 220<^^> 

5.52 

2200 

-18.2 

4.12 

700 

-18.0 

2^ 

-2.8356 ± 0.0003^“^ 

0.0 

-2.84 

- 

-12.6 

-2.85 

- 

-12.1 

2^ 

2.76±0.1l(“> 

800 ± 130(“> 

1.75/2.36 

375/1150 

-12.6/-3.0 

1.72/2.35 

161/920 

-12.1/-3.0 


2.86 ± 

2080 ± 300('^1 







2.t 

3.88 ±0.05^“^ 

430 ± 80^“> 

3.90/6.87 

163/2010 

-12.6/-3.0 

5.04/5.48 

570/3600 

-12.1/-3.0 


6.0±0.2('=1 

1700 ± 200^'=^ 

5.33 

3700 

16.3 

5.48 

3800 

13.4 


6.81 ± 0.02l“l 

258 ± 15^“> 

6.82 

620 

16.3 

6.81 

800 

13.4 

43 + 



13.1 

1800 

16.3 

13.1 

2060 

13.4 


(a) Ref. [3; (b) Ref.^; (c) Ref.[||; (d) Ref.d], (e) Ref. [Hi 


width of 0.258 MeV has been known for a long time 
(Ref.0). Much more recently, an additional 4+ state 
has been reported with an energy of 6.0 MeV (above 
threshold) and a width of 1.7 MeV [l^. Our calcula¬ 
tions are consistent with the existence of a well-defined 
and relatively narrow resonance, which by means of the 
three-body force can be made to appear at 6.81 MeV, in 
agreement with the experimental value for the 4^ state. 
The computed width for this resonance is 0.62 MeV with 
the Ali-Bodmer potential and 0.80 MeV with the Buck 
potential. When this is done it can be seen that a quite 
broad resonance (^3.7 MeV wide) appears numerically 
with both potentials at and energy of about ^5.4 MeV. 
This broad resonance can be interpreted as the one given 
experimentally in Ref. M- Also, both the Ali-Bodmer 
and Buck potentials, give rise to a dg state at 13.1 MeV, 
having also a similar width in the vicinity of 2.0 MeV. 

As seen in Table im the energies and widths of the 
resonances are rather independent of the choice made 
for the potential. In general, with the two potentials 
used, a similar three-body force gives rise to pretty much 
the same resonance properties. However, there are a few 
exceptions, mainly the 0;^ and 2'^ states, and the narrow 
2^ resonance at 1.7 MeV, which with the Ali-Bodmer 
potential is more than a factor of two broader than with 
the Buck potential. These differences can lead to visible 
differences in the differential cross sections, and therefore 
also in the transition strengths, for the cases where these 
potential-dependent resonances are involved. 

As an illustration of how the resonance properties have 
been obtained, we show now in Fig.[T]the discretized con¬ 
tinuum spectra after a complex scaling calculation for the 
O’*" (Fig.[TJi), 2+ (Fig.IlJ)), and 4+ (Fig.[T}:) states in 
We show only the results obtained with the Buck poten¬ 
tial. The calculations have been done using a complex 
scaling angle, 6, of 0.30 rads for the 0+ and 4+ states. 


and 0.25 rads for the 2+ states. In all the three cases 
the energy cuts starting from the origin, from the res¬ 
onance in ®Be, and from the 4+ resonance in ®Be are 
clearly seen. Note that, due to the very small energy of 
the 0+ resonance in ®Be, the energy cut starting from the 
0 "^ state overlaps with the cut starting from the origin, 
which corresponds to strict three-body continuum states. 

These cuts are rotated in the complex energy plane by 
an angle 29, and the states in the cuts correspond to pure 
continuum states. Typically, the states corresponding to 
the three-body resonances fall clearly out of the energy 
cuts, and their position is independent of the scaling an¬ 
gle 9. These resonances are indicated in the figure with 
the corresponding labels (the bound 0)*' and 2^ states 
are not shown in the hgure). However, there are sev¬ 
eral cases where the resonances lie pretty close to the 
continuum states, and their identification is not so obvi¬ 
ous. The natural way to isolate these resonances from the 
continuum states would be to increase the value of the 
scaling angle 9. However this is not always an efficient 
solution, due to the technical difficulties arising from the 
use of a too large scaling angle. For narrow resonances 
close to the threshold, like the Hoyle state (0^ state), 
the resonance wave function falls exponentially to zero 
very fast, and the resonance can easily be identified by 
looking directly into the resonance wave function. For 
not very narrow resonances, like the 2^ state in Fig. [Hd, 
or the 4^ in Fig.[TJ;, it is much more efficient to separate 
the resonance from the continuum states by increasing 
the attraction of the three-body force, and trace how the 
resonance moves when the attraction of the three-body 
force is progressively released to the desired value. This is 
illustrated in Figs. [TJ) and Figs. [TJ:, where the red crosses 
indicate the position of the resonances when the strength 
of the three-body force is put equal to —12.1 MeV in the 
2+ case and 7 MeV in the 4+ case. The red arrows show 
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TABLE III. Computed values of prms = (in fm) and r^ms (in fm) for the same states shown in Table HIl 



prms 

f^rms 

prms (Buck) 

Tr-ma (Buck) 

Ot 

6.9 

2.5 

6.8 

2.5 

02+ 

11.5 

3.6 

11.4 

3.6 

O.t 

11.3-t 14.0 

3.5-t 11.1 

10.9-^12.3 

3.4-k 10.6 

2t 

6.8 

2.5 

6.7 

2.4 


9.8-t 13.9/10.7-t 12.4 

3.2-h 11.0/3.6 -t 10.7 

9.9-^12.9/10.0 -t 12.6 

3.2-t 10.7/3.2-t 10.7 

2^ 

6.0-t 13.1/6.9-^10.9 

2.2-h 10.7/2.5 -t 10.2 

8.3-h 11.0/7.5 -h 16.8 

2.8-t 10.2/2.4-t 11.7 

4t 

9.3 -h 11.7 

3.1 -h 10.4 

9.5 -t 11.6 

3.1-k 10.4 

4 

6.0-h 10.1 

2.3-t 10.1 

5.9-t 10.1 

2.3-k 10.1 

4 

7.7 -h 11.6 

2.7-h 10.4 

7.9 -t 11.5 

2.7-k 10.4 


how the resonances move when when decreasing the at¬ 
traction. Using this procedure the 2'^ and 4f resonances 
can be unambiguously identified. 

It is interesting to note that the states denoted as 4^ 
and 4^ in Fig. [T]: do actually cross when decreasing the 
repulsion of the three-body force (red crosses in the fig¬ 
ure), in such a way that eventually, for a sufficiently large 
three-body attraction, the narrow resonance becomes the 
first 4+ state and the broad one becomes the second 4+ 
state. In other words, when decreasing the repulsion of 
the three-body force, the state denoted in Fig. [TJ; as 4j 
appears closer to the threshold than the 4f state. It is 
important to keep this fact in mind in order to determine 
which of the two first 4+ states should be assigned to the 
first band, and which one to the second. 

C. C radii 

In order to get a feeling of the size of the states, we 
show in Table Ell the expectation value for the 

same cases shown in Table HIl This expectation value will 
be denoted as Prms, and it is computed within the com¬ 
plex scaling frame, meaning that p^ms = ('I'|p^e®^®|'I'), 
where 41 is the complex rotated wave function of the res¬ 
onance and 6 is the complex scaling angle used in the 
calculation m- Therefore, the expectation value is not 
necessarily real even though the p-coordinate is. With 
this value it is possible to obtain the root mean square 
radius of each state {r^ms in Table Ell), which is given 
by: 

rrms = + Rl, ( 18 ) 

where Ra = 1.47 fm is the root mean square radius of 
the a-particle. Again, due to the complex scaling cal¬ 
culation, the computed rrms of a resonance is in general 
a complex number. As discussed in Ref.jl^, the imag¬ 
inary part of a computed observable (the energy is the 
most obvious example) can be attempted interpreted as 
the uncertainty of the value given by the real part. This 


quantity permits a fast comparison between the spatial 
extensions of the different states. 

As seen in the table, the results obtained with the Ali- 
Bodmer and Buck potentials are very similar. It is re¬ 
markable that, on the one side, the 0)'' and 2f states have 
similar sizes, both of them in the vicinity of Prms ^ 6.8 
fm and rrms ^ 2.5 fm, and, on the other side, the same 
happens with the 0^ and 2^, resonances, with a value of 
Prms 10 fm and rrms ^ 3.4 fm. This may suggest that 
these two sets of 0+ and 2+ states could correspond to 
two different rotational bands each of them with a rea¬ 
sonably well-“frozen” structure. It is interesting to note 
that for the 4^ and 4^ states the computed values of 
Prms are precisely about 10 fm and 6 fm, respectively, 
which might indicate a crossing of the first and second 
4+ states, in such a way that the 4+ state could belong 
to the second band and the 4^ state could belong to the 
first one. This is consistent with the discussion in Fig.E;, 
where we showed that the 4^ state becomes actually the 
first 4+ and the 4^ state becomes the second 4+ when 
the repulsion in the three-body force is diminished. 

Finally, a third band could be present containing the 
0;^, 2;^, and dg resonances. However, the values of prms 
given in Table ED for these three states is not very stable 
(it ranges from 6.0 fm to 11.3 fm), and furthermore, as 
seen in Table El the experimental energy of the 0;^ state 
is higher than the one of the 2j^ state, which makes it 
very unlikely that these two states belong to the same 
rotational band. 


V. ELECTRIC QUADRUPOLE CROSS 
SECTIONS 

In this section we shall consider electric quadrupole 
cross sections for transitions between the different O'*", 
2+, and 4+ states in These transitions involve bound 
states and resonances. As already mentioned, in this 
work we do not treat the resonances separately. We just 
compute the wave functions for the discrete continuum 
states obtained by imposing a box boundary condition. 
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FIG. 1. (Color online) Fourth quadrant of the energy plane 
showing the discrete O’*', 2^, and 4'*' continuum spectra in 
after a complex scaling calculation with the Buck potential. 
The angle 9 on each panel is the scaling angle used in the 
corresponding calculation. The strength of the three-body 
force is —18.0 MeV, —3.0 MeV, and 13.4 MeV for the 0’^, 
2^, and states, respectively. The bound 0)'' and 2^ states 
are not shown. The lowest resonances for each of the three 
angular momenta are indicated by the corresponding labels. 
In the 2'^ and 4'^ cases the red crosses indicate the position 
of the resonances when the strength of the three-body force is 
reduced from —3.0 MeV to —12.1 MeV, and from 13.4 MeV to 
7 MeV, respectively, and the arrows show how the resonances 
move when decreasing the three-body attraction. 

The only role of the resonance energies obtained in the 
previous section is to provide information about where 
to put the energy windows when computing the cross 
sections. 

In general, the cross section, as a function of the in¬ 
cident energy, is given by Eq.®, where the summation 


over j is restricted to the states within the chosen final 
energy window. For clean isolated resonance peaks it 
is easy to choose a meaningful, and rather well defined, 
window around the peak. For more complicated energy 
dependent cross sections this could easily be more am¬ 
biguous. However, this would only reflect the more com¬ 
plicated physical nature resulting from overlapping and 
perhaps interfering resonances. Variation of choices of 
windows and subsequent analyses are then related to re¬ 
duction of contributions to properties of individual reso¬ 
nances. Thus the ambiguity in the choice of the window 
still remains but now containing information about the 
underlying physics 

We shall here define the windows as ± T/j, where 
Eh and T h are the energy and the width of the resonance 
in the initial or final state (we shall use the computed 
energies and widths given in Table El. Obviously, when 
the final state is a bound state the corresponding width is 
zero, and the summation over j in Eq.® disappears and 
the wave function in the right part of the matrix element 
refers to the bound-state wave function. 

The continuum states have been discretized by impos¬ 
ing a box boundary condition to the radial solutions fnip) 
in the set of Eqs. m- In particular, a box size of 200 fm 
has been used, which amounts to an average energy sepa¬ 
ration of about 0.03 MeV between the lower-lying contin¬ 
uum states. This means that for an energy window with 
a width in the scale of MeV we will have a significant 
amount of discrete continuum states within the final en¬ 
ergy window. However, this procedure can not be applied 
to transitions into the Hoyle state. This resonance is ex¬ 
tremely narrow (about 8 eV), which implies that only a 
huge box could provide a significant number of discrete 
continuum states within a final energy window of only a 
few eV wide. 

Discretization of the continuum by use of such a huge 
box of perhaps 30000 fm compared to the 200 fm used 
in this work is clearly meaningless within the present nu¬ 
merical approach. The reasons are the huge number of 
necessary partial waves, the increasing basis size for each 
of the corresponding components, and the coupling of 
all the potentials due to the long-range Coulomb force. 
Therefore the transitions into the Hoyle state will be 
treated as ordinary transitions into a bound state. On 
the other hand this is a very accurate approximation for 
such a small width. 


A. 2^ —0’*' transitions 

In Fig 12] we show the cross section for the a + a-\- a ^ 
-|- 7 reaction for a transition between the continuum 
2+ states and the ground state in The three-body 
force for the O'*" state has been chosen to reproduce the 
experimental separation energy of into three alphas 
(see Table nil). The figure shows the results obtained with 
the Ali-Bodmer potential (solid curves) and the Buck 
potential (dashed curves). Two different (structureless) 
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transition in the a + Q + a—^ ^^C + 7 reaction as a function of 
the three-body energy. The results with the Ali-Bodmer po¬ 
tential (solid cnrves) and the Buck potential (dashed curves) 
are shown. The (a) and (b) panels correspond to the calcu¬ 
lations placing the 2"^ resonance at 1.7 MeV and 2.4 MeV, 
respectively. 


three-body potentials have been used to place the 2 + 
resonance at different energies. 

In panel (a) the 2+ energy is chosen to reproduce the 
measured binding energy of the bound 2^ state. This 
implies that the first resonance, 2 ^, is found at about 
1.7 MeV, which gives rise to the pronounced peak in the 
cross section around this energy. As we can see, the cal¬ 
culation using the Buck potential produces a taller peak 
than with the Ali-Bodmer potential. This is related to 
the fact that the computed 2 ^ state is clearly narrower 
when the Buck potential is used. By fitting these peaks 
with the Breit-Wigner shape in Eg. lfTOl) we extract the 
value of r.y for the 2 ^ —>■ O)*" reaction, which for the 
Ali-Bodmer and Buck potentials is found to be 205 meV 
and 115 meV, respectively. We emphasize that, as seen 
from Eq. (USD, in the center of the Breit-Wigner shape 
{E = Eb) the value of the cross section is proportional 
to r.y/r_R (assuming that T.^ <C E/j), and therefore dif¬ 
ferent values of the cross section in the peak does not 
necessarily imply different values of the r.y-width, or, in 
other words, only for similar values of F/j the difference in 
the peaks of the cross section is directly translated into a 
difference in the T-y-values. The computed r.y-widths are 
collected in Table El Note that the presence of the 2^ 
resonance at about 4 MeV with the Ali-Bodmer potential 
or 5 MeV with the Buck potential has no visible effect 
on the cross section. In fact, the small peak observed 
at around 6 MeV is basically produced by the broad 2^ 
resonance, which with both potentials appears close to 6 
MeV, as shown for the Buck potential by the red crosses 
in Fig. [Ud. 

In the second calculation (Figl^Ja) the three-body force 
for the 2 + states has been chosen such that the lowest 2 + 


TABLE IV. T-^-widths, in meV, for the resonance-to- 
resonance transitions obtained after fitting the peaks in the 
cross sections for the different reactions with the Breit-Wigner 
shape given in Eg. (1161) . For the reactions where the 2^ res¬ 
onance enters, the two resonance energies, 1.7 MeV and 2.4 
MeV, have been considered. 


-£■2+ ^ 

Ali-Bodmer 

1.7 MeV 2.4 MeV 

Buck 

1.7 MeV 2.4 MeV 

0 

t 

-flM 

205 


160 

115 

175 

2+ 0+ 

1.0 


5.9 

0.4 

4.5 

2+^0+ 

— 


1950 

— 

4300 

23+ ^ 0+ 

0.7 


80 

6.5 

190 

2+ -)■ 2+ 

6.4 


22 

4.0 

22 

2+^2+ 

20 


320 

88 

950 

4+^2+ 

130 


44 

118 

57 

4j 2+ 

36 


25 

26 

18 

4+^2+ 


50 


50 

4+ ^ 2+ 


635 


860 

4+^2+ 


3100 


3500 


resonance is placed at about 2.4 MeV, in closer agree¬ 
ment to the experimental value given in Ref.@. The 
peaks in the cross sections are therefore shifted towards 
higher energies. As we can see, the new peaks are now 
clearly broader and almost a factor of 10 lower than be¬ 
fore. Again, the calculation with the Buck potential gives 
rise to a taller peak than when the Ali-Bodmer potential 
is used. The r-,,-widths obtained with the two poten¬ 
tials are in this case very similar, 160 meV and 175 meV, 
respectively. In this case the second peak in the cross 
section located around 7 MeV is produced by the 2)^ 
resonance. From this peak we can also extract the 
widths for the 2 ;^ —>■ 0 )/' reaction, which for completeness 
are also given in Table HVl 

In Ref.Q the F.y-width of the 2^ resonance in has 
been extracted from the photo-dissociation cross section 
of the O)/' ground state, and it has been found to be 60 ± 
10 meV, which is roughly a factor of 3 smaller than the 
values obtained in this work when the 2^ resonance is 
placed at 2.4 MeV (which is the energy that agrees better 
with the one given in the same reference). In Ref.@, 
where the same experiment is analyzed, the F-^-width of 
the 2j state is quoted to be 135 ± 14 meV, which agrees 
better with the results of our computation. However, 
as mentioned above, the relevant quantity when fitting 
the data with a Breit-Wigner shape is F.y/F/j, which in 
both, Ref.0 and Ref.@, is of about 6.5 • 10“®. This 
value is roughly a factor of 3 smaller than the value we 
have obtained with the 2^ resonance at about 2.4 MeV 
(F-y/Fi{ « 17 • 10“®). Therefore the photo-dissociation 
cross section for the ^^C-l -7 —>■ a + a + a shown in 
Refs.jl, 1^ should also be roughly a factor of 3 smaller 
than the one obtained in this work. 

The photo-dissociation cross section can be extracted 
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FIG. 3. (Color online) Photo-dissociation cross section for the 
^^C-l -7 a-|-a-|-a reaction between the 0 ^ ground state and 

the continuum 2 "*" states as a function of the photon energy. 
The experimental curve has been multi^ied by a factor of 
3. The experimental data are from Refs.Q (squares) and 
(circles). Only the results with the 2"^ resonance at 2.4 MeV 
are shown for both, the Ali-Bodmer (solid) and Buck (dashed) 
potentials. 


from the cross sections shown in Figj^] simply by using 
the relation in Eq.(j4]). The result is shown in Figl3l where 
the experimental data (multiplied by a factor of 3) are 
the ones given in Refs.Q (squares) and 0 (circles) and 
the curves are our calculations with the Buck and Ali- 
Bodmer potentials (and the energy of the 2^ resonance 
at about 2.4 MeV). As we can see, the shape of the exper¬ 
imental cross section is well reproduced. The disagree¬ 
ment on the absolute value of the cross section may be 
related to the fact that only a fraction of the 0 + ground 
state in corresponds to the three-alpha structure used 
in this work. Another possibility is that improvements of 
the experimental analyses with the available conflicting 
results [1,0 also would change the normalization. 

In FigUlwe show the same as in Figl2]for transitions 
from the 2+ continuum states in into the Hoyle state 
(treated as a bound state). Again, the (a) and (b) panels 
show the result when the 2^ resonance is at 1.7 MeV and 
2.4 MeV, respectively. The computed cross sections are 
clearly smaller than for transitions into the 0 ^ ground 
state. Together with the main peak produced by the 2^ 
resonance, several other peaks produced by higher 2 + 
states are clearly seen. In particular, in FigH^ a small 
peak can be seen at an energy of around 4 MeV for the 
Ali-Bodmer potential and 5 MeV for the Buck potential. 
These two peaks are produced by the 2^ resonance given 
in Table im An additional peak produced by the 2'^ res¬ 
onance at about 6 MeV is also seen in the cross section. 
The computed F-^-widths obtained after fitting the lower 



FIG. 4. (Golor online) The same as Fig[2] for the 2^ —^ Oj 
transition in the a + a + a ^ ^^G-l -7 reaction. 

peak in FigH^ with the Breit-Wigner shape in Eq. (HU) 
are 1.0 meV and 0.4 meV with the Ali-Bodmer and Buck 
potentials, respectively. 

When the 2^ resonance is at 2.4 MeV (FiglUa) the 2'^ 
state becomes very broad and moves up in energy, giv¬ 
ing rise to the corresponding bumps in the cross section 
observed in EigHjD. In the case of the Buck potential a 
bump produced by the 24 state at about 11 MeV (FiglTJo) 
is also seen. The computed F.^ widths for the transition 
from the 2^ state to the Hoyle state are in this case 5.9 
meV and 4.5 meV with the Ali-Bodmer and the Buck 
potentials, respectively. For completeness, we also give 
in Table IS the computed values of F-^ for the 2 ;^ —>■ 0 ^ 
transition. 


B. 4+ —2"*’ transitions 

In FigO we show the cross section for the 4+ —>■ 2f 
transition. The three-body force in the 2+ states is such 
that the binding energy of the 2^ bound state matches 
the experimental value. The two potentials, Ali-Bodmer 
and Buck, give very similar results with a very sharp peak 
at 6.8 MeV, which is actually produced by the 4j reso¬ 
nance. The effect of the broad 4^ state can be seen in the 
little shoulder that appears in the cross section around 
5 MeV. This is more clearly seen in FiglHl where we fit 
for the Ali-Bodmer (panel (a)) and Buck (panel (b)) po¬ 
tentials the cross section peak as a linear combination of 
two Breit-Wigner functions with the parameters for the 
4^ (dashed curves) and 4^ (dotted curves) resonances 
given in Table UTl The sum of the two Breit-Wigner curves 
is shown by the thin-solid curves, which match pretty 
well the peak of the computed cross sections (thick-solid 
curves). The coefficients in this linear combination give 
directly the F^y widths for the 4f —>• 2f and 4^ —)> 2f 
transitions, which are, with the Ali-Bodmer and Buck 
potentials, 50 meV and 80 meV for the first transition. 
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FIG. 5. (Color online) The same as FiglJ] for the —>• 2f 

transition in the a + a + a —>■ + 7 reaction. 



FIG. 7. (Color online) The same as Fig[2] for the 4'*' —^ 2^ 
transition in the a + o + a—>■ ^^C + 7 reaction. In panel (b) 
the thin-solid and thin-dashed curves are the cross sections 
for the Ali-Bodmer and Buck potentials, respectively, using a 
cutoff for the photon energy of 2 MeV. 



E (MeV) E (MeV) 


FIG. 6. (Color online) Fit of the sharp peak in the cross sec¬ 
tion for the 4'^ —>■ 2f transition by means of two Breit-Wigner 
functions with the parameters corresponding to the 4^ and 
4^ resonances given in Table HIl Panels (a) and (b) show the 
fit when the Ali-Bodmer and Buck potentials, respectively, 
are used. The thick-solid curve curve is the computed cross 
section as shown in FigO The dashed and dotted curves are 
the Breit-Wigner functions corresponding to the 4f and 4^ 
resonances, respectively, whose sum is shown by the thin-solid 
curve. 


and 635 meV and 855 meV for the second transition, re¬ 
spectively. The broad peak observed at about 13 MeV 
is mainly produced by the 4^ resonance, which has very 
similar properties with the two potentials. The T..^ values 
for the 4^ —>■ transitions are also given in Table HVl 

In Figl3 we show the cross sections for the 4+ —>■ 
transitions. Again, two possible energies, 1.7 MeV and 
2.4 MeV, have been considered for the 2^ resonance (pan¬ 
els (a) and (b), respectively). This is a transition into a 
resonant state, which implies that the summation over j 



E (MeV) E (MeV) 

FIG. 8. (Golor online) Fit of the cross section for the 4^ —^ 2^ 
transition by means of two Breit-Wigner functions with the 
parameters corresponding to the 4^ and 4^ resonances given 
in Table im Panels (a) and (b) show the fit when the Ali- 
Bodmer and Buck potentials, respectively, are used and the 
2+ resonance is found at 1.7 MeV. Panels (c) and (d) show 
the same as (a) and (b) when the 2^^ resonance energy is 
2.4 MeV. The thick-solid curve is the computed cross section 
as shown in Fig[7] (with the photon energy cutoff in panels 
(c) and (d)). The dashed and dotted curves are the Breit- 
Wigner functions corresponding to the 4f and 4^ resonances, 
respectively, whose sum is shown by the thin-solid curve. 

in Eq. (jS]) is made over the discrete continuum final states 
within the energy window E/jiT/j, where and T/j are 
the energy and width of the final resonance. Therefore, 
the value of the cross section will be fully determined 
by the arbitrary choice for the width of the final energy 
window. 

The calculation of the continuum-to-continuum cross 
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section presents also the numerical complication of the 
so-called infrared catastrophe (see Ref.jl^), which ap¬ 
pears in the region of small values {E close to E') 
due to the l/E^ dependence of the bremsstrahlung cross 
section at small photon energies. Therefore, the cross 
sections shown in Figl7]by the thick-solid (Ali-Bodmer) 
and thick-dashed (Buck) curves will be contaminated by 
this effect, also called soft-photon contribution, for ini¬ 
tial energies {E) in the vicinity of the final energy win¬ 
dow (centered at ^ 1.75 MeV in panel (a) and at ^ 2.4 
MeV in panel (b)). A fully relativistic treatment of the 
bremsstrahlung cross section would correct this anomaly 




In Fig[7^ the effect of the soft-photon contribution is 
seen in the little shoulder observed at about 2 MeV. This 
peak is therefore unphysical and it does not correspond to 
any resonance in the initial state. The other two peaks, 
at about 5.5 MeV and 7.2 MeV, correspond to the 4^ and 
4^ states, although the large width of the 4f resonance 
gives rise to a large interference between the two lowest 
4+ states. 


In the case shown in Figl7)3 (2^ resonance at ^ 2.4 
MeV) the final energy window reaches a value of around 
3.5 MeV. For this reason the soft-photon peak and the 
one corresponding to the broad 4f resonance can not be 
distinguished in this case. Usually the soft-photon contri¬ 
bution is removed by introducing a cutoff in the photon 
energy. In FiglT^ the thin-solid and thin-dashed curves 
are the cross sections with the Ali-Bodmer and Buck po¬ 
tentials, respectively, when a cutoff of 2 MeV is used for 
Ej. This cutoff has been chosen in order to place the 
4^ peak at an energy of ^ 5.5 MeV, similar to the one 
observed in panel (a), that we know from the value of 
the 4+ energy should be the correct position. Although 
the choice of the cutoff value is to some extent arbitrary, 
it should not be much bigger than 2 MeV, since the en¬ 
ergy separation between the 4^ and the 2j resonances in 
FigEla is of about 3 MeV. For this reason, the transition 
strengths obtained from the lowest peak in Fig[7}D (thin 
curves) for the 4f —>• 2j process should be taken as a 
minimum value. A choice of the cutoff energy smaller 
than 2 MeV would increase the strength. The second 
peak at ^ 7.1 MeV is produced by the 4^ state. The 
position and height of this peak is not affected by the 
cutoff in the photon energy, but, again, it contains an 
important interference from the 4f resonance. 

In order to extract the r.y-values for the 4'^ —^ 2^ and 
4j —>■ 2^ transitions, and due to the large interference 
between the two lowest 4'^ resonances, it is convenient to 
proceed as discussed in FiglHl and make a simultaneous 
fit of the 4^ and 4^ peaks in the cross sections by means 
of Eg. lfTOl) . These hts are shown in Fig[51 The panels (a) 
and (b) are the calculations with the Ali-Bodmer and 
Buck potentials for the 2+ state at 1.7 MeV, and panels 
(c) and (d) are the same calculations with the 2+ state at 
2.4 MeV. In these two last panels the ht has been made 
to the cross sections obtained after removing the soft- 
photon contribution by means of the cutoff in E^ (thin 




FIG. 9. (Color online) The same as Fig[2] for the 2^ —^ 
transition in the a + a + a ^ ^^C-l-y reaction. 


curves in Fig[7}D). The Breit-Wigner curves correspond¬ 
ing to the d)*" and 4^ states are shown by the dashed and 
dotted curves, respectively, and they have been obtained 
using F-y-values equal to 130 meV and 36 meV, respec¬ 
tively, in panel (a), 118 meV and 26 meV in panel (b), 
44 meV and 25 meV in panel (c), and 57 meV and 18 
meV in panel (d). The sum of the two Breit-Wigner fits 
is given by the thin-solid curves, which clearly disagree 
with the computed total cross section (thick-solid curves) 
for large energies. This is due the fact neither the contri¬ 
bution from the wide 4^ state given in Table |lT] nor the 
contribution form the continuum background have been 
included in fit. 

Finally, we show in Figj^the corresponding cross sec¬ 
tions for the 2+ —>■ 2^ transitions. The peaks produced 
by the 2^ and 2^ resonances are clearly seen in all the 
cases. The F-^-widths for the 2^ —>■ 2^ and 2j^ —^ 2f 
transitions with the two potentials and the two different 
energies of the 2^ resonance considered in this work are 
again given in Table HVl 


VI q ( e 2 ) transition strengths 

The photo-emission processes from a bound state into 
other bound states or into the continuum is described un¬ 
ambiguously by the transition strength. This is not true 
any more for continuum to continuum transitions even 
when rather well-defined resonance peaks are present 
in the corresponding cross sections. These transition 
strengths are decisive indicators for structure similari¬ 
ties between excited states, where prominent examples 
are collective rotational or vibrational states built on the 
same intrinsic configuration. In our present case of three 
alpha-particles such collective rotations have been sug¬ 
gested many times. This necessarily involves continuum 
structures with all the related ambiguities. Still, we want 
to know if collective rotations is a reasonable description 
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TABLE V. B{E2) transition strengths in for different transitions between states. For transitions from continuum 

states, B~i and B^ denote the strengths obtained from Eq. 113 and from the total (integrated) cross section, respectively. 
For the reactions involving the 2^ resonance the two cases, with the resonance energy at 1.7 MeV and 2.4 MeV, have been 
considered. The second, third, and fourth columns give the known experimental values, and the results obtained with the 
microscopic cluster model calculation described in Ref.[l3| (labeled as a-cluster) and with the Antisymmetrized Molecular 
Dynamics (AMD) method. The corresponding references are given in the table. 


Transition 

Exp. 

a-cluster 

AMD 

Ali-Bodmer 

Buck 

2+^0+ 

7.6 ±0.4^“) 

9.16^“^ 

8.4(d) 


10.2 




9.9 


4+^2+ 



15.8(‘^> 

Bj = 9.0 

Ba 

= 8.2 

By 

= 12.5 

Ba 

= 9.9 





A2+=1.7 MeV 

E2+ = 

2.4 MeV 

E 2 + = 

1.7 MeV 

E 2 + = 

2.4 MeV 





B^ Ba 

By 

Ba 

By 

Ba 

By 

Ba 

2+^0+ 



102(d) 

230 175 

170 

140 

120 

135 

158 

136 

4+ ^ 2+ 



595('^) 

170 135 

153 

130 

150 

130 

220 

145 

2+ ^ 0+ 

2.6±0.4i“) 

0.84i“) 

5.l(‘*( 


0.91 




0.96 


4+^2+ 




B-y = 1.2 

Ba 

= 0.6 

By 

= 1.2 

Ba 

= 1.0 





A2+=1.7 MeV 

E2+ = 

2.4 MeV 

E2+ = 

1.7 MeV 

E 2 + = 

2.4 MeV 





By Ba 

By 

Ba 

By 

Ba 

By 

Ba 

2+ 0+ 

0.73 ±0.13'*'^ 
1.57±0.13(‘^> 

1.99^“^ 

0.4(d) 

4.1 3.1 

2.0 

1.9 

2.4 

2.8 

2.3 

2.2 

2+ —> 2^ 

^2 ^ 




3.6 2.9 

5.5 

5.2 

2.5 

2.9 

5.5 

5.4 

4t 2+ 



7.5(‘*( 

8.5 4.0 

9.8 

5.0 

6.5 

3.5 

7.5 

4.7 


(a) Ref.0; (6) Ref.jl]; (c) Ref.i]; (d) Ref.[il 


of some of the states in the ^^C-spectrum. We therefore 
first discuss the detailed numerical results from our three- 
body model, and in the following subsection we relate to 
the simplest possible rotational model. 


A. Numerical results 

In Sect. IIII Cl we described the two methods used in 
this work to extract the transition strength for reactions 
involving continuum states. In the first method we inte¬ 
grate the cross section (Eq.Q) over two energy windows 
chosen around the initial and final resonance energies. 
In this work we have taken the windows as En ± F^j 
where F/j is the width of the resonance. This procedure is 
equivalent to integration of the cross sections computed 
in the previous section over the initial energy E under 
the peaks corresponding to the initial resonant state. 
The total cross section obtained in this way is divided 
by the constants multiplying the transition strength (see 
Eqs.® and d?]))- These transition strengths will be de¬ 
noted by Ba-- The application of this method requires 
well-separated peaks in the cross sections for each of the 
resonances, in such a way that the integration under a 
given peak contains very little contamination from the 
interference with other states. However, this is not al¬ 
ways the case in our calculations. In particular, the 4^ 
and 4^ states are very close to each other, giving rise 
to a large interference between them, as seen in EigslB] 
and m For this reason, for transitions having the 4^ 


or 4^ resonances as initial state, the integration of the 
cross section over the initial energy window will be made 
considering not the computed cross section (thick solid 
curves in Figs E] and |8|), but the Breit-Wigner curves cor¬ 
responding to each of the resonances (dashed curves for 
the 4^ state and dotted curve for the 4^ state). In the 
second method we make use of Eq. (HU, where F.^ 
been extracted after fitting the computed cross section 
with the Breit-Wigner shape in Eg. dTSl) . The values of 
F.^ have been given in Table ITVl The transitions strengths 
computed in this way will be denoted by 


The first two lines in Table in show the computed 
strengths for the 2f —>• 0)^ and 4^ —>• 2^ transitions, 
which are tempting to associate with transitions between 
states belonging to the first rotational band in In 
here we have taken into account that, as suggested by the 
computed values of p^ms (see Table iBll) and the behavior 
of the 4^ and 4^ resonances when modifying the three- 
body force (FiglT]), the 4^ state is expected to belong to 
the first band and the 4f state to the second. The first 
transition in the table involves only bound states, and 
the results obtained with the Ali-Bodmer and Buck po¬ 
tentials are similar, 10.2 e^fm^ and 9.9 e^fm'^ which are 
in good agreement with the result in Ref.[T^ from a mi¬ 
croscopic a-cluster model calculation. These values are a 
bit bigger than the one obtained in Ref. 1^ where the mi¬ 
croscopic Antisymmetrized Molecular Dynamics (AMD) 
method was used. The results are also bigger than the 
experimental value [T5 | which most temptingly can be at¬ 
tributed to a better match between the 3a wave functions 
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than nature prescribes. 

For the 4^ —>■ transition the computed ,B.y-strengths 

agree slightly better with the AMD calculation than the 
Scr-values, which are a bit smaller. In comparison to 
other model calculations it is important to know the con¬ 
ceptual difference between our continuum calculation and 
these bound-state treatments. When a resonance is lo¬ 
cated in the continuum with a substantial width it is 
very likely that the transition strength is reduced due 
to a smaller overlap of wave functions than by assuming 
one bound-state like resonance wave function. It is also 
important to keep in mind that the choice of the energy 
windows when computing is arbitrary, and a small 
increase in such windows will give rise to values of 
closer to Bj. 

The next two lines in the table give the i?2-transition 
strengths between the states in the second sequence of 
energies, i.e., the 2^ —>■ 0^ and 4:^ —^ 2^ transitions. In 
order to see the dependence on the energy of the 2 J reso¬ 
nance, we show, as in Table ITVl the results corresponding 
to a 2^ energy of 1.7 MeV and 2.4 MeV. 

For the 2j —>■ 0^ transition, when the 2^ resonance 
is at 1.7 MeV, the results with the Buck potential are 
clearly smaller than the values obtained with the Ali- 
Bodmer potential. This is particularly true for B-y, where 
the difference is of almost a factor of 2. This is due to the 
fact that the width of the 2^ resonance (at 1.7 MeV) with 
the Buck potential is around half the width of the one 
with the Ali-Bodmer potential (see Table ID). As seen in 
FigHh-, the height of the cross section peak corresponding 
to this resonance is similar for both potentials. This im¬ 
plies that r.y/rfl is also similar in both cases. Therefore, 
a Ffl-value about a factor of 2 smaller gives rise to a F-^- 
value also about a factor of 2 smaller, and consequently, 
as seen in Eq. a transition strength also about a fac¬ 
tor of 2 smaller. For a 2^ energy of 2.4 MeV, which is 
in better agreement with the experimental value given 
in [1, 0, the computed strengths with the Ali-Bodmer 
and Buck potentials are now closer to each other (the 
difference between the corresponding F/j values is now 
small). Also, the By and Ba values are reasonably con¬ 
sistent with each other, and they are clearly bigger than 
the AMD result (102 e^fm"^) given in Ref. [l^. In any 
case this difference is relatively unimportant compared 
to the sensitivity of the computed transition strengths 
on the methods used. 

For the if —>■ 2^ transition the results obtained with 
the two potentials are consistent with each other, even 
if By changes from 153 e^fm'* to 220 e^fm^ when the 2^ 
state is at 2.4 MeV. One has to take into account that, 
together with the inherent uncertainties in the determi¬ 
nation of By and B^, for reactions involving the 4^ and 
4^ states we have to deal with the additional uncertainty 
arising from the interference between the two resonances 
and the soft-photon contribution. For this reason, the 
agreement between the By and B^ values seen in Ta¬ 
ble |V] for the 4^ —>■ 2J reaction can be considered quite 
acceptable. In any case, all the results given in the Table 


for this reaction are significantly smaller than the result 
given in Ref. 0. 

In the lower part of Table El we show the computed 
transition strengths between states belonging to differ¬ 
ent bands. For the —>■ 0^ transition (where the Hoyle 
state is treated as a bound state) we obtain a strength 
of ~ 0.9 e^fm'^, in reasonably good agreement with the 
result of the microscopic a-cluster calculation given in 
Ref. [T^, but clearly smaller than the experimental value. 
For the 2^ —)> O)*" transition our computed strengths 
are quite stable, although the agreement with the mi¬ 
croscopic a-cluster calculation given in Ref.fl^ is bet¬ 
ter, as expected, when the energy of the 2^ resonance 
is put at about 2.4 MeV, also in better agreement with 
the experimental value. However, these energies of about 
2 MeV are in clear disagreement with the AMD result. 
When compared to the experimental data, our strength is 
about a factor of 3 bigger than the result given in Ref. ^ , 
but in better agreement with the value given in Ref. 9|, 
where the same experiment as in Ref.[^ is reexamined, 
see FigH For the 4^ —>■ 2^ transition our results are rea¬ 
sonably consistent with the AMD calculation, especially 
when considering the By values. The computed strengths 
for the transitions 4f —>■ 2^ and 2^ —^ 2^ are also given, 
although for these cases the comparison with previous 
results or experimental data is not possible. 

In general, the computed Ba- and By strengths are con¬ 
sistent with each other, especially if we take into account 
that the computed Ba values are obviously dependent 
on the width chosen for the energy windows, and the 
By strengths are very sensitive to the value of Ey used 
in Eq. ClI, since the photon energy appears to the fifth 
power. This is particularly true when the initial and fi¬ 
nal energies are not very far from each other, like for 
instance in the 2^ —>■ Oj transition with the 2^ energy 
at 1.7 MeV. In this case a variation of Ey from 1.3 MeV 
to 1.4 MeV (and F-^ = 1 meV) implies a change in By 
from 230 e^fm^ to 330 e^fm^. Together with this, the re¬ 
sults involving the 4^ and 4^ states have the additional 
uncertainty arising from the interference between them 
and the soft-photon contribution. 

When comparing the results with the Ali-Bodmer and 
Buck potentials, the general conclusion is that there are 
no significant differences between them. Only in the tran¬ 
sitions involving the 2^ resonance at 1.7 MeV an impor¬ 
tant difference, mainly for the By values, has been ob¬ 
served. This is due to the very different width obtained 
for this resonance with each of the potentials. Therefore, 
from the transition strengths it is not possible to answer 
the question of what potential is more appropriate in or¬ 
der to describe the alpha-alpha interaction. 


B. Rotational model 

For transitions between states within a schematic ro¬ 
tational band, and assuming axial symmetry for the sys- 
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tern, the quadrupole transition strength is given by (4ll | 
^ J') = y|^Oo(^0; 20| J'0)2, (19) 

where the projection, K, of the angular momentum on 
the intrinsic symmetry axis has been assumed to be zero. 
The intrinsic quadrupole moment Qo is given by: 

Qo = -xf- yj)) , (20) 

i 

where i runs over all the charged particles with charge 
Qi and whose center of mass coordinates are denoted by 
{xi,yi,Zi), where the 0 -axis is chosen along the intrinsic 
symmetry axis. The expectation value is taken in the 
intrinsic body-fixed coordinate system. The quadrupole 
moment is a measure of the deformation and it has ideally 
one characteristic value for a sequence of states belonging 
to a given rotational band. From each of the calculated 
transition strengths given in Table El it is then easy, to 
obtain from Eq. m the absolute value of the intrinsic 
quadrupole moments. The |(5o|-values obtained in this 
way are given in Table ED 

The intrinsic quadrupole moment Qo is related to the 
static quadrupole moment Q of each individual state with 
angular momentum J by the simple expression: 

Q = ^J^{JJ\d,o\JJ) = ( 21 ) 

where the operator 0\^ is given in Eq.® and, again, a 
band with K = 0 has been assumed. 

Eor the transitions between the states in the first band, 
2'^ —Oj*" and 4^ —>• 2'^ (the state 4^ is the one assigned 
to the first band), the computed values of IQol in TableEJ 
are rather stable, consistent with each other, and consis¬ 
tent as well with the experimental value and previous 
calculations. The value of the static quadrupole moment 
Q can be easily computed for the bound 2^ state, and 
we have obtained 6.6 efm^ and 6.5 efm^ with the Ali- 
Bodmer and Buck potentials, respectively. These values 
agree with the experimental value of 6 ± 3 efm^ given 
in [ 4 ^. From Eq. m we can then extract the intrin¬ 
sic quadrupole moment for the 2^ bound state, which is 
found to be —22.7 efm^ and —23.0 efm^ with the Ali- 
Bodmer and Buck potentials, respectively. This result 
agrees as well with the —21.6 efm^ given in the Ref. [4^ . 
These values are also consistent with the ones quoted 
in the first two rows in Table Ell for the 2f —>• Oj*" and 
4j —>• 2f transitions. These results support the idea of 
the states in this band as arising from the rotation of a 
given intrinsic state, and also the fact the 4^ state is the 
one actually belonging to the first band. 

Eor the transitions in the second band, 2^ —>■ 0^ and 
4)*" —>■ 2^, the relative difference between the computed 
values of |Qo| is higher than the one found in the first 
band. Typically, the results for the 4f —2^ transition 
are even a 25% smaller than the ones for the 2^ —>■ Oj 
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FIG. 10. (Color online) Coordinates of the three-alpha system 
in the center of mass frame. 


transition. If we restrict ourselves to the results involv¬ 
ing the 2^ state at 2.4 MeV, which is in better agree¬ 
ment with the experimental value, the computed |Qo| 
values range from ^ 70 efm^ up to ^ 90 efm^. Taking 
into account all the uncertainties already discussed, espe¬ 
cially when the 4f state is involved, we can consider that 
these results are consistent with each other. In this con¬ 
nection, it is important to remember that, as discussed 
in FiglTlo, the transition strengths obtained from the 4f 
peak for the 4'^ —>■ 2^ reaction should be taken as its 
minimum value. A small decrease in the photon energy 
cutoff would enhance the 4f peak in EiglTlD, making the 
strengths of the 4f —>• 2^ process closer to the ones of 
the 2^ —Oj reaction. Therefore, we can conclude that 
the transition strengths given in Table I VII between the 
states in the second band in are also consistent with 
the idea of a relatively well ’frozen’ structure and the ro¬ 
tational character of the band. For the same reason, we 
can also consider that our results are not dramatically 
far from the 72 efm^ and 145 efm^ obtained in Ref. (T^ . 

In the lower part of Table I VII we give the intrinsic 
quadrupole moments for the transitions between states 
in different bands. It is interesting to note that for the 
4'^ —>■ 2^ transition, which in principle should be a tran¬ 
sition belonging to the first rotational band, we obtain a 
value for |Qo| even a factor of 3 smaller than the values 
quoted in the upper part of the table for the 2f —>■ 0)'' 
and 4^ —>■ 2'^ transitions. This result is then also consis¬ 
tent with the assignment of the 4j state to the first band 
and the 4^ state to the second. In fact, as also seen in 
the table, the |Qo| value for the 4^ —>• 2^ is also clearly 
smaller than the values for the transitions between the 
states in the second band. 

A different approach can be made by applying the def¬ 
inition given in Eq. (1201) to a system made of three point¬ 
like Qf-particles in a schematic triangular structure, as 
shown in EigfTUl From the coordinates given in the fig¬ 
ure, we can obtain the hyperradius, psb, of such a three- 
body system, which is given by: 

Pi. = ^y(3 + 4(d/a)^), 


( 22 ) 
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TABLE VI. Absolute value of the intrinsic transition quadrupole moments IQol (iu efm^) obtained from the B{E2) transition 
strengths given in Table IVl and Ea. (ll9ll . \Qo\cr and IQ 0 I 7 denote the intrinsic quadrupole moments obtained from the transition 
strengths Ba and B-y, respectively. 


Transition 

Exp. 

a-cluster 

AMD 

Ali-Bodmer 

Buck 

2 +^ 0 + 

19.5 ±0.5 

21.46 

20.5 



22.6 




22.3 


4+^2+ 



23.6 

IQoIt = 17-8 

IQol 

= 17.0 

IQol 

7 = 21.0 

Qo c 

= 18.7 





^2+ ~ 

1.7 MeV 

-^ 2 + ~ 

=2.4 MeV 

£^ 2 + = 

=1.7 MeV 

^2+ ~ 

2.4 MeV 





Qo 7 

\Qo\t7 

IQ0I7 

IQolo- 

IQ0I7 

IQolo- 

Qo 7 

|Qo|o 

2 + -^ 0 + 



72 

108 

94 

92 

84 

78 

82 

89 

83 

4+ ^ 2+ 



145 

77 

69 

73 

68 

73 

68 

88 

71 

2 +^ 0 + 

11.4 ±0.9 

6.50 

16.0 



6.8 




6.9 


4+^2+ 




IQol 

7 = 6.5 

Qo 

=4.6 

Qo 

7 = 6.5 

IQol 

0 = 5.9 





E2+ = 

1.7 MeV 

E2+ = 

=2.4 MeV 

E 2 + = 

=1.7 MeV 

^2+ ~ 

2.4 MeV 





Qo 7 

\Qo\t7 

IQ0I7 

IQolo- 

IQ0I7 

|Qo|o 

Qo 7 

|Qo|o 

2 + 0 + 

6.03 ±0.54 

10.0 

4.5 

14.3 

12.5 

10.0 

9.8 

11.0 

11.9 

10.8 

10.5 


8.88 ±0.37 











2 +^ 2 + 




11.2 

10.1 

13.9 

13.5 

9.4 

10.1 

13.9 

13.8 

4 ^ 2 + 



16.2 

17.3 

11.9 

18.6 

13.3 

15.1 

11.1 

16.2 

12.9 


where ttIq is the mass of the a-particle and m is the 
normalization mass used to define the Jacobi coordinates 
(the nucleon mass in our calculations). 

Also, using the coordinates given in FiglTOl we can get 
the moments of inertia relative to each of the three co¬ 
ordinate axes (where axis-3 is perpendicular to the plane 
containing the three particles shown in the figure). These 
three moments of inertia take the form: 

Ti = (23) 

12 = 2ma(f (24) 

3 

1 3 =Ii +I 2 — -maO? + 2ma(f, (25) 

where in li the index i refers the axis with respect to 
which the moment of inertia has been calculated. 

We have to keep in mind that in Eqs. m and o 
an axial symmetry has been assumed. This implies that 
these equations can be used in those cases in which two 
of the moment of inertia in Eas. ((23|) to dSSl) are equal 
to each other, while the axis with respect to which the 
moment of inertia is different to the other two defines the 
symmetry axis (and therefore the z-axis in Eg. (1201) I. In 
other words, there are three possible geometries for the 
system in Fig ITU] to which Eas. dTUl) and (lUU)) can be ap¬ 
plied. They are the geometries corresponding to having 
the symmetry axis along axis-1, axis-2, or axis-3: 

(a) Symmetry axis along axis-1 (I2 = ^ 3 )- Making 
equal Eas. dMl) and (jUUl) we get that this happens when 
a = 0, leading to X 2 =X^ = 2mad? and Xi = 0. Making 
a = 0 corresponds to the geometry shown in FigfTTh. 
where the three alphas are aligned along axis-1, which is 
the symmetry axis. Therefore when computing Qo by use 
of Eg. ljUUl) the z-axis has to be taken along axis-1, and 


axis-2 

(a) 

I2=l3 

±1=0 


axis-1 

symmetry 

axis 



, , axis-2 

(c) 

, A =±2 


axis-1 

,• 

/:^axis-3 

1 symmetry 
. 1 axis 


FIG. 11. (Color online) Geometries in the three-alpha system 
in Fig. 1101 corresponding to the symmetry axis along axis-1 
(o), along axis-2 ( 6 ), and along axis-3 (c). The third axis, 
axis-3, is perpendicular to the plane holding the three alphas. 


for the three-alpha system {qi = 2 e) we get Qo = Sed^, 
which by means of Eo. dUUl) and taking rria/m k, 4 can be 
written as Qo = epgj,. 

(b) Symmetry axis along axis-2 {Xi = X 3 ). Making 
now equal Eas. dUUl) and (IU51) we get that this happens 
when d = 0, and we get in this case Ii = I3 = 3771^0^/2 
and I2 = 0. The corresponding geometry is now the one 
shown in FigfTTb. with the three particles along axis-2, 
one of them at the position a and the other two at —a/2. 
Taking then axis-2 as the z-axis, we get from Eg. dUUl) 
that for the three alphas Qo = 6ea^, which making use 
of Ea. (l22|) can again be written as Qo = epg^. 

(c) Symmetry axis along axis-3 (Ii = X 2 ). As be¬ 
fore, making equal Eos. (l23l) and (l24t we find that the 
symmetry axis is along axis-3 when (Pjo? = 3/4, and 
Xi = X 2 = 3TOQa^/2, and I3 = 3maa^. In this case 
the corresponding geometry is the one shown in Fig lllb . 
which corresponds to the three alphas in an equilateral 
triangle. Taking then the z-axis along axis-3 (perpen¬ 
dicular to the plane holding the three alphas) we get 
from Eg. (|UU)) that Qq = —6ea^, or, using again Eo. dUUl) . 
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Qo — ^p\b/‘^- 

Summarizing, the three axially symmetric geometries 
shown in FiglTTllead to either Qo = which happens 
for the linear geometries (a) and (&), or Qo = 
which happens for the equilateral triangular geometry 
(c). 

As already discussed, our three-body calculations are 
consistent with Qo ~ —23 efm^ for the bound state, 
as given for instance in [d^, |4^ . Such a negative value 
for the intrinsic quadrupole moment is only consistent 
with Qo = —ep 3 f,/ 2 , and therefore with the equilateral 
structure. Furthermore, as shown in Table Ell for the 
states in the first band we have that Prms ~ 6.7 fm, from 
which we get Qo = Ri —22.5 efm^. This result 

is in very good agreement with the computed intrinsic 
quadrupole moment for the state, the values given in 
[^1^, and also with the values given in Table IVll for the 
transitions between the states in the first band. We then 
conclude that the computed transition strengths and in¬ 
trinsic quadrupole moments for the transitions between 
the states in the first band in the three-alpha system are 
consistent with the ones corresponding to transitions be¬ 
tween the states of a rotational band in a three-alpha 
system where the three alphas are sitting in the vertices 
of an equilateral triangle. 

For the states in the second band, we have from Ta¬ 
ble nni that Prms ~ 10 fm. Using this value we get 
that for the linear structures in FigfTTk and FigfTTb 
the intrinsic quadrupole moment in the band should be 
Qo ~ 100 efm^, and for an equal sided triangular distri¬ 
bution (FigfTTbl it should be Qo ~ —50 efm^. Keeping 
in mind that, as discussed in connection with FiglTjo, 
the results shown in Table EH for the 4:f —^ 2 ^ transi¬ 
tion could have been underestimated, we can conclude 
that our computed intrinsic quadrupole moments for the 
states in the second band are consistent with the ro¬ 
tational estimate of ^ 100 efm^ corresponding to the 
aligned structure shown in either FigfTTk or FigfTTb. 

Another hint about the rotational character of the two 
bands in can be given by the sequence of energies 
in the { 0 f, 2^,42 } and the {Oj, 22 ,4^} bands. As it 
is well known, in the case of axial symmetry and K = 0 
bands, these energies should follow the rule: 

Ej-Eo = ^J{J + l), (26) 

where Ej is the energy of the state in the band with 
angular momentum J, Eq is the energy of the lowest 
state in the band, and I is the moment of inertia relative 
to an axis perpendicular to the symmetry axis. 

For the three axially symmetric configurations given in 
Fig HD we have: 


Geometry (a): 

X = 2maCp = maplbl4: 

(27) 

Geometry ( 6 ): 

X — ‘imatP 12 = maP3j/4 

(28) 

Geometry (c): 

X = 3maa^/2 = rfiap\bl'i, 

(29) 

where, again, we have made use of Eg. (1221) with a 

= 0 in 


case (a), d = 0 in case ( 6 ), and cP/a^ = 3/4 in case (c). 


For the first band, for which Eq = —7.28 MeV, the 
experimental energies of the 2 f and 4^ states given in 
Table El lead by means of Ea. (l26l) to fP/2X k, 0.74 MeV 
and 0.71 MeV, respectively. These values are very similar 
to each other, supporting the fact that the first band 
fulfills the condition of a “frozen” structure that gives 
rise to a rotational band. Furthermore, this first band has 
been seen to be consistent with the equal sided triangular 
structure and Prms ~ 6.7 fm. Making use of Ea. (l2^ 
we estimate from the rotational model that fpI2X k, 0.9 
MeV, which is relatively close to the values obtained from 
the experimental 2 /" and 4/" energies. 

Eor the second band, for which Eq = 0.38 MeV, the 
energies of the 2 J and 4/" states in Table El determine 
fPI2X R! 0.33 MeV and 0.25 MeV for the 2^ —>■ 0^ and 
4/" —>■ 2^ transitions. From the analysis of the intrin¬ 
sic quadrupole moments we concluded that the states in 
the second band could at best be consistent with the 
linear structure in Figslllb and lllb . which by use of 
Eg. dTfl) or (1^ . and taking prms ~ 10 fm (as shown in 
Table ED, lead to the estimate from the rotational model 
fPI2X R 0.2 MeV, which is again reasonably consistent 
with the values obtained from the energy differences. 
Thus, also in this case the analysis of the energies in the 
{Oj, 2^,4/"} band supports the conclusion obtained from 
the intrinsic quadrupole moments, namely, the properties 
of the states in the second band of the three-alpha sys¬ 
tem are consistent with the behaviour expected for states 
belonging to a rotational band whose structure would be 
the one in FigfTTk or FigfTTb. 

Although the analysis of both, Qo and I, does not per¬ 
mit to distinguish between the linear structures shown in 
FigsfTTk andfTTb. it is important to note that the alpha- 
alpha Coulomb repulsion clearly prevents the formation 
of a system as the one described in Figlllb. For this 
reason, together with the fact that the geometry given in 
FigfTTb is not consistent with cluster model calculations, 
we can conclude that the only possible aligned structure 
for the states in the second band in is the one given 
in FigfTTk. 

We emphasize that our rotational model analysis is 
schematic, even within our a-particle model. The O'*" and 
2 ^ bound states are known only to be approximately de¬ 
scribed by an a-particle model. The continuum states 
investigated in the present paper are on the other hand 
expected to be better a-particle states. However, this 
does not imply that a schematic classical linear or trian¬ 
gular structure provide an accurate description. Already 
the quantum mechanical probability distributions differ 
from this simplified picture. Both an equal sided triangle 
and a linear chain can only reveal the qualitative essence 
of the structure of these states. In particular, the ac¬ 
curacy of this illuminating qualitative picture is not in 
conflict with the bent arm structure of the Hoyle state 
found in many previous calculations as well as in the 
present one. 

The present interpretation of a rotational band struc¬ 
ture is at first glance in conflict with the interpretations 
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of both the O t a nd 2^ states as a very dilute gas of three 
particles [3 ^|- These two interpretations are difficult 
to reconcile, since a gas of particles does not have ro¬ 
tational states as defined in classical textbooks [s^ and 
used in the present paper. Also, the radii of these states 
in [3,14^ differ apparently from each other, although as¬ 
cribed to different energies but with “similar structure”. 
In any case, our conclusions about energies, transition 
probabilities, and quadrupole moments rely entirely on 
the computed wave functions. In this sense, interpreta¬ 
tions do not play any role, since they do not enter into 
any of the quantum mechanical solutions. Whether these 
solutions are understood in terms of a frozen intrinsic 
structure or a dilute gas-like structure is not important, 
and both interpretations can in principle be correct si¬ 
multaneously. 


VII. SUMMARY AND CONCLUSIONS 

In this work we investigate the spectrum of and 
its rotational character. The existence of two O'*", 2'^, 
sequences, related to the 0^ ground state and the 0^ res¬ 
onance (the Hoyle state), led people to refer to these two 
series as rotational bands in The sequence, 0+, 2+, 
4+, suggests that the geometry corresponds to axial and 
i ?2 symmetric intrinsic states. However, transitions be¬ 
tween rotational states build on different intrinsic shapes 
with K = 0 bandheads would also be mathematically 
described in precisely the same way. This means that 
the same tests apply even if other states with different 
quantum numbers are members of the same band. 

In the present case we investigated whether the prop¬ 
erties of these continuum states are consistent with the 
electromagnetic transition strengths for a rotational band 
built on the ground and first excited O'*" states. In partic¬ 
ular, the transition strength should equal the square of 
the intrinsic quadrupole moment multiplied by a Clebsch- 
Gordan coefficient. 

The three-body wave functions in have been ob¬ 
tained through the hyperspherical adiabatic expansion 
method, and the continuum spectrum has been dis¬ 
cretized by imposing a box boundary condition on the hy¬ 
perradius. The discretized continuum spectrum has been 
obtained on the real energy axis, without preferential 
treatment of the resonances. Thus, we have from the out¬ 
set the natural distribution over continuum states. Most 
existing other models treat resonances as bound sates 
even though many of them have a substantial width. All 
transition strengths are then present in one wave bound 
state like resonance function. 

The transition strength for the different reactions has 
been extracted from the corresponding y-emission cross 
sections. For transitions into a resonance, the calculation 
of the cross section as a function of the initial three-body 
energy requires specification of the range of energies in¬ 
vestigated in the final state, which should be around the 
final resonance energy. For this reason, some informa¬ 


tion about the position of the resonances is needed. In 
this work this is done by means of the complex scaling 
method, which provides the resonances as poles of the 
iS-matrix. It is important to emphasize that this is done 
only to know what energy windows to consider for tran¬ 
sitions into a resonance, since the method used to dis¬ 
cretize the continuum and to compute the cross sections 
only deals with discretized continuum states with real 
energies. 

From the cross sections two different methods have 
been used to obtain the transition strengths. The first 
one uses the F.y-width as a parameter to fit the peaks of 
the cross section, associated to resonances in the initial 
state, with the usual Breit-Wigner shape. From the F.^- 
width the transition strength is computed making use of 
the simple expression relating these two quantities. In 
the second method the cross section, as a function of 
the incident energy, is integrated under the peak corre¬ 
sponding to initial resonance. This integrated cross sec¬ 
tion gives, except for some known factors, the transition 
strength. 

The cross sections for the 4+ to 2+, the 2+ to O'*', 
and 2+ to 2+ transitions have then been computed. Two 
different potentials, the Ali-Bodmer and the Buck po¬ 
tentials, have been used to describe the a-a interaction. 
The main features of the computed cross sections are in 
general independent of the potential used. Only those re¬ 
actions involving resonances whose properties are more 
sensitive to the potential show a sizable difference in the 
cross section. This happens for instance in the reactions 
with initial states with angular momentum and parity 
2 +, especially when the lowest 2+ resonance is located 
at about 1.7 MeV, for which the predicted resonance 
width differs by a factor of 2 depending on the potential 
used. In any case, when the energy of the 2^ resonance 
is moved up to 2.4 MeV (in better agreement with the re¬ 
cent experimental value) the dependence on the potential 
is much less relevant. 

The transition strengths obtained from the cross sec¬ 
tions are in general consistent with each other. First, the 
results are similar no matter what method is used to ex¬ 
tract the transition strengths or especially when 
taking into account the uncertainties inherent to each of 
the methods. And, second, they are also independent of 
the potential used. Again, only in the reactions involving 
the 2^ state at 1.7 MeV show a higher dependence on the 
potential. 

The first result we obtained is that the 4f and 4^ 
states, which belong to different bands in have ac¬ 
tually crossed, in such a way that the 4)*" state belongs 
to the second band and the 4^ to the first one. This 
is first suggested by the fact that the computed values 
of Prmsi which should be similar for systems having the 
same ’frozen’ spatial structure, provides values for the 
4^ state similar to the ones of the Oj and 2^ states, and 
values for the 4^ state similar to the ones of the 0)'' and 
2^ states. Furthermore, investigating how the 4f and 
4^ resonances move in the complex energy plane when 



21 


making the effective three-body force more and more at¬ 
tractive, we have seen that the first 4+ state becoming 
bound would actually be the 4^. 

The assignment of the 4f state to the second band 
and the 4^ to the first is also confirmed by the analysis 
of the transitions strengths and the intrinsic quadrupole 
moments. The values of Qq are obtained assuming a 
rotational character for the two bands in under in¬ 
vestigation, in such a way that the transition strength 
is basically the square of the intrinsic quadrupole mo¬ 
ment multiplied by some geometrical factor depending 
on the initial and final angular momenta of the transi¬ 
tion. Doing like this we have seen that for the 4^ —>■ 2^ 
reaction the computed transition strength is clearly more 
consistent with the previous AMD calculation than the 
one corresponding to the —>■2^. Furthermore, the in¬ 
trinsic quadrupole moments derived from the transition 
strengths for the —>■ 0^ and 4j —>■ reactions are 

rather similar and consistent with previous calculations 
and the experimental value. This consistency would dis¬ 
appear if the reaction 4f —?> 2 '^ were the one taken into 
account. The same happens for the reactions in the sec¬ 
ond band. Although the computed transition strength 
for the 4^ —>■ 2^ is even a factor of 3 smaller than the 
previous AMD calculation, this discrepancy is certainly 
less important than when considering the 4^ —>■ 2 ^ re¬ 
action, and it can actually be understood from the un¬ 
certainties associated to the method used to extract the 
strength. Also, the intrinsic quadrupole moments for the 
2 j ^ Oj and 4f —7> 2^ are reasonably stable, stability 
that would be clearly broken if assuming the 4j —>■ 2j 
transition as belonging to the second band. Therefore, 
from the stability of the intrinsic quadrupole moments 
for each of the {0^ ,2^ ,4^} and {0j,2^,4^} bands we 
can conclude that each of the bands correspond to states 
having a rather well preserved rigid structure. 

The rotational character of two bands is confirmed 
when comparing with the prediction obtained from an ax¬ 
ially symmetric rotating structure made by three point¬ 
like alpha particles. Assuming values for p^b similar to 
the ones computed numerically for the states in the first 
and second bands {prms ^ 6.7 fm and ^ 10 fm, respec¬ 
tively) we have seen that the Qo values previously com¬ 
puted for the states in the first band are consistent with a 
equilateral triangular structure rotating around an axis 
perpendicular to the plane holding the three particles. 
For the states in the second band, the computed Qo val¬ 
ues are mostly consistent with a linear distribution of 
the particles. The same results are obtained from the 
analyses of the moments of inertia. 

Summing up, we have calculated genuine continuum- 
continuum electric quadrupole {E2) transitions from 0'*', 
2+ and 4+ states in 12C. The transition strengths are 
not well defined since neither initial nor final states have 
well defined bound-state like structures in the contin¬ 
uum. Energy windows around resonance positions ex¬ 
tending at least the natural widths must contribute to the 
transitions. Instead we design, approximate and deduce 


correspondingly non-observable transition strengths from 
(double) differential cross sections. The results are two 
sequences of rotational band structures with ground and 
Hoyle states as band heads with interchange of the order 
of the two close-lying 4+-resonances. The derived intrin¬ 
sic quadrupole moments and moments of inertia are con¬ 
sistent with axial triangle and almost linear structures, 
respectively. In conclusion, our procedure to compute 
continuum-continuum transitions are used to classify six 
12C-states in two rotational bands. 


Appendix A: Three-body incoming flux 

Given the three-body reaction a-|-6-l-c —?> A-I- 7 , 
the cross section corresponding to this process requires a 
definition of the incoming flux of particles. This can be 
done similarly to the two-body case, where the incoming 
flux is well defined. 

At the two body level the Schrodinger equation takes 
the form: 

84/ n 

where p, is the reduced mass, 4' the two-body wave func¬ 
tion, V,. is the three-dimensional gradient operator in 
terms of the relative vector r, and V is the two-body 
potential. As described in any Quantum Mechanics text 
book, from the Schrodinger equation we can obtain the 
flux of incoming particles as: 

j = ^ - (Vr4'*)4'], (A2) 

2 pi 

which for an incoming plane wave leads to the well-known 
result 

j = (A3) 

with p being the two-body relative momentum. 

At the three-body level, the Schrodinger reads: 

where now px is the reduced mass of two of the parti¬ 
cles, is the gradient operator in terms of the relative 
coordinate Tx between those two particles, py is the re¬ 
duced mass of the third particle and the two-body system 
made by the other two, and ^ry is the gradient operator 
in terms of the relative coordinate Vy between the third 
particle and the center of mass of the other two. The 
wave function 4' is now a three-body wave function, and 
V contains all the potentials involved in the three-body 
system. 

Let us int roduce now t he usu al Jacobi coordinates 
X = yjPx/tTiXx and y = ^Py/mry which contain 
the arbitrary normalization mass m. From this defini¬ 
tion we immediately get that = \/Px/rn^x and 
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'Vvy = y^/iy/mVy, where 'Vx and Vy are the gradient 
operators in terms of the x and y Jacobi coordinates, 
respectively, and which permit to write the Schrddinger 
equation Ea. (IA4l) as 


.^ 5 ^ T 

'"-at = 


V'i>, 


(A5) 


which is equivalent to Eq. m and where we have intro¬ 
duced V2 ^ = V2 + V 2 

Therefore, from Ea. dASI) . and proceeding exactly as 
done at the two-body level, we can define the three-body 
incoming flux as: 


J ^ ['I'* - (V.,yVE-*)vI/] , 


(A6) 


where ^x.y is the six-dimensional gradient operator in 
the Jacobi coordinates x and y. 

Considering now the incoming plane wave 
we get the incoming flux: 


flux of outgoing particles normalized with the incoming 
flux, and therefore, according to Eq. (lAlOl) , the cross sec¬ 
tion will depend as well on m. In other words, the cross 
section is only well defined for a given definition of the 
hyperradius. 


Appendix B: Cross section for the a + b + c^A + 'y 
process 

As discussed above, the cross section for the process 
a + b + c ^ A-I -7 depends on the normalization mass m. 
It is not a well defined physical observable. In fact, the 
usual observable (therefore independent of m) for this 
kind of reactions is the reaction rate, which after multi¬ 
plication by the density of particles gives the number of 
reactions per unit time and unit volume in some specific 
environment. The reaction rate for the a + b + c ^ A-I -7 
process at a given three-body kinetic energy E, Rabc{E), 
is given in Eq.(l) of Ref.js^: 


j = (A7) 

m 

where {k^, ky) represents the six-dimensional momentum 
given by the three-body momenta and fey, which are 
the momenta associated to the Jac obi coordinates x and 
y. If we define now n = we then get: 

j = ft-, (A8) 

m 

which is completely analogous to the two-body result in 
Eg.dMt. 

The result given above is correct provided that the 
wave function is normalized to 1 in terms of the Jacobi 
coordinates x and y. However, the correct normalization 
should be done not in terms of the Jacobi coordinates, 
but in terms of the relative coordinates and Vy. Since 

J d£cdy|4'p = ^ ^ J dvxdryl'i’l'^, (A9) 


we have that to get normalized to 1 in terms of the 
relative coordinates we have to multiply by the factor 
{mj that when done in Eq. (IMl) leads 
to the final expression for the incoming flux: 




(AlO) 


It is important to note that the flux of particles given 
above depends on the arbitrary normalization mass m. 
This is reflecting the fact that for a given value of the hy¬ 
perradius p = the relative distances between 

the three incoming particles, and Vy, will be different 
for different values of the normalization mass. It is then 
clear than the flux of particles through some given hy¬ 
persurface has to depend as well on the choice of m. The 
cross section for the reaction a-|-6-l-c—^-A-l-yis just the 


Rabc{E) = (Bl) 

, ft3 Stt 2{2Ja + 1) (E-y 

^ {2Ja + l){2Jb + 1)(2J, + 1) \ ¥ 

where v is the number of identical particles in the three- 
body system, J^, Jf,, and Jc are the angular momenta of 
particles a, ft, and c, Ey is the photon energy, and cfy{Ey) 
is the photo-dissociation cross section of the process A + 

7 —>■ a J- ft J- c. 

The reaction rate is given by the product of the cross 
section and the incoming flux. Using the expression in 
Eg. dAlOp we can then write: 

i?, 6 ,(E;)=ft-(—) ^ aabc{E), (B2) 

m \p,xj KPyJ 

which after use of Eq. m leads to the following relation 
between the cross section for the process a-|-ft-l-c—>■ A-I -7 
and the one corresponding to the inverse reaction: 

cJabcjE) ^ ^ , 2(2Ja + 1) 32^ 

ay{Ey) {2Ja + l){2Jb + l){2J^ + l) \nc) ’ 

(B3) 

which depends on the normalization mass m through k. 


^ ay{Ey), 


Appendix C: Reduced matrix element 


In this work the three-body wave functions are written 
as given in Eq. ®: 

1 


4'.7 = 


,5/2 




(Cl) 


The angular wave functions 0;((p, H) in Eo. dCip are 
expanded in terms of the hyperspherical harmonics: 


n) = Yy c^^Hp) ® x: 


s 

Ssc S'u 


(C2) 
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where q collects the quantum numbers 
{K,£x,£y,L,Sx,S}, where K is the hypermomen¬ 
tum, £x and Sx are the relative angular momentum and 
total spin of the two particles connected by the x Jacobi 
coordinate, £y is the relative momentum between the 
third particle and the center of mass of the first two, L is 
the total angular momentum obtained by coupling of £x 
and £y, and S is the total spin obtained by coupling Sx 
and the spin of the third particle Sy. The total angular 
momentum J is obtained by coupling of L and S. The 
total spin function of the three-body system is denoted 
by Xs^sy^ the hyperspherical harmonics yxe^iyi^) 
are defined as: 

x[Yi^{nx)(8)Yi^{nyX^X (C3) 


where v is such that K = 2v + £x + £y and the normal¬ 
ization constant can be found for instance in [^ . 

For obvious reasons, in order to compute the matrix 
element ||ypTA(cp)||‘hj/^) it is convenient to write 
the three-body wave function in terms of the Jacobi set 
{a;p,yp}. When this is done, insertion of Eas. dCjp and 
(IC2I) into Ea. (jCll) gives the full expansion of the three- 
body wave function. Using this expression the integra¬ 
tion over rix and fly involved in the calculation of the 

matrix element {XXlVpYXXl^j^M') can be trivially 
made, and only the integrals over p and a remain. 

After some algebra, we get the final result: 




X£y£'yLL'JJ' 


L X L' 

£' £ £ 


q n' q' 

-X, 


J J' A [ (-1) 

V L S \ 2“+''+2 


\/4^ 


£y X £'y\ 

0 0 0 y 


(C4) 


where the indeces without and with primes refer to the 
<i)j^ and wave functions, respectively, u = y/2u Y 1, 
a = £x Y 1/2, and b — {X + 1 + £y + £'y)l2.). Finally, 

^(xiyt'KK' , 1-1 ,1 , 

denotes the integral over the hyperangle a. 


which can be written as: 


(.xiyl'yKK' 

X 


(C5) 



dx{l - a;)“(l + 
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